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Abstract. The theoretical and experimental results concerning the thermodynamical and low-frequency trans- 
port properties of hybrid structures, consisting of spatially-separated conventional low-temperature supercon- 
ductor (S) and ferromagnet (F), is reviewed. Since the superconducting and ferromagnetic parts are assumed 
to be electrically insulated, no proximity effect is present and thus the interaction between both subsystems is 
through their respective magnetic stray fields. Depending on the temperature range and the value of the exter- 
nal field Hext, different behavior of such S/F hybrids is anticipated. Rather close to the superconducting phase 
transition line, when the superconducting state is only weakly developed, the magnetization of the ferromagnet 
is solely determined by the magnetic history of the system and it is not influenced by the field generated by the 
supercurrents. In contrast to that, the nonuniform magnetic field pattern, induced by the ferromagnet, strongly 
affect the nucleation of superconductivity leading to an exotic dependence of the critical temperature Tc on Hext- 
Deeper in the superconducting state the effect of the screening currents cannot be neglected anymore. In this 
region of the phase diagram various aspects of the interaction between vortices and magnetic inhomogeneities 
are discussed. In the last section we briefiy summarize the physics of S/F hybrids when the magnetization of 
the ferromagnet is no longer fixed but can change under the influence of the superconducting currents. As a 
consequence, the superconductor and ferromagnet become truly coupled and the equilibrium conflguration of 
this "soft" S/F hybrids requires rearrangements of both, superconducting and ferromagnetic characteristics, as 
compared with "hard" S/F structures. 

Some figures in this paper are in color only in the electronic version. 
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List of main notations 

Acronyms 

GL Ginzburg-Landau, 

DWS domain-wall superconductivity, 

F fcrromagnet or ferromagnetic, 

OP order parameter, 

RDS reverse-domain superconductivity, 

S superconductor or superconducting 

ID onc-dimensional 

2D two-dimensional 

Latin letters 

A vector potential, corresponding to the total magnetic field: B = rot A, 

a vector potential, describing to the nonuniform component of the magnetic field, b = rot a, 

B total magnetic field: B = Hext + b, 

b nonuniform component of the magnetic field induced by ferromagnet, 

c speed of light, 

Ds thickness of the superconducting film, 

Df thickness of the ferromagnetic film (or single crystal), 

/ absolute value of the normalized OP wave function: / = y^(Re'^)^ -|- (Im^)^, 

jext the density of the external current: rotHe^t = (47r/c) je^t, 

js the density of superconducting currents, 

Gsf free (Gibbs) energy of the S/F hybrid, 

Gm term in the free energy functional accounting for the spatial variation of the magnetization, 

Hext external magnetic field. 

Hex exchange field, 

Hci lower critical field: H^i = $o ln(A/^)/(47rA2), 

Hc2 upper critical field: Hc2 = $o/(27r^^) = H^°^ (1 - T/Tco), 

H^2^ upper critical field at T = 0: H^^^ = $0/(277^2)^ 

h separation between superconducting and ferromagnetic films, 

L angular momentum of Cooper pairs (vorticity): = f{r) e^^"^ , 

e.H magnetic length: = ^/W^J(2tt\H^^ , 

il effective magnetic length determined by a local magnetic field 6*: il = y/ $o/(27r|6j|), 

ijp typical width of the localized OP wave function, 

M magnetization of the ferromagnet, 

Ms magnetization of the ferromagnet in saturation, 

mo dipolar moment of a point- like magnetic particle, 

Rg radius of the superconducting disk, 

Rf radius of the ferromagnetic disk-shaped dots, 

Rrf position of a point-magnetic dipole: = {Xj^, Y^, Z^}, 

Tco superconducting critical temperature at 5 = 0, 

w period of the one- dimensional domain structures in ferromagnet 

Greek letters 

a, f3 constants of the standard expansion of the density of the free energy with respect to |^'|^, 

el^-* sclf-cncrgy of the vortex line per unit length: = {^q/AttX)'^ InA/^, 

8 the OP phase: O = arctan(Im^/Re^/'), 

A temperature-dependent magnetic field (London) penetration length: A = Ag / a/1 — T/T^q, 
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Ao magnetic field penetration length at T = 0, 

^ temperature-dependent superconducting coherence length: ^ = ^o/v^l — T/Tcq, 

^0 Ginzburg-Landau coherence length at T = 0, 

TT 3.141592653..., 

p electrical resistivity, 

$0 magnetic flux quantum: <I>o = irhc/e ~ 2.07 Oe-cm^, 

superconducting order parameter (OP) wave function, 

OP saturated value, "^q = ^/—a/|3, 
ij) normalized OP wave function, V = ^/^o 

Coordinate systems 

Throughout this paper we use both cartesian reference system (x, y, z) and cylindrical reference system 
(r, ip, z), where z— axis is always taken perpendicular to the superconducting film/disk. 
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1. Introduction 

According to the classical Bardeen-Cooper-Schrieffer theory of superconductivity, the ground state of the 
superconducting condensate consists of electron pairs with opposite spins (the so-called spin-singlet state) 
bounded via phonon interactions [UIl]- As early as 1956, Ginzburg [3] pointed out that this fragile state of 
matter could be destroyed by the formation of a homogeneous ferromagnetic ordering of spins if its corresponding 
magnetic field exceeds the thermodynamical critical field of the superconductor. Later on, Matthias et al. 
[H [51 [S] demonstrated that besides the orbital effect (i.e. a pure electromagnetic interaction between the 
ferromagnetic and superconducting subsystems), there is also a strong suppression of superconductivity arising 
from the exchange interaction which tends to align the spins of the electrons in detriment of Cooper-pair 
formation. Anderson and Suhl [7] predicted that a compromise between these antagonistic states can be achieved 
if the ferromagnetic phase is allowed to break into domains of size much smaller than the superconducting 
coherence length ^ in such a way that from the superconductivity point of view, the net magnetic moment 
averages to zero. Alternatively, Larkin and Ovchinnikov [8j and Fulde and Ferrel [9], theoretically predicted 
that superconductivity can survive in a uniform ferromagnetic state if the superconducting order parameter is 
spatially modulated. 

In general terms, the effective polarization of the conduction electrons, either due to the external field 
Hext (orbital effect) or the exchange field Hex (paramagnetic effect), leads to a modification (suppression and 
modulation) of the superconducting order parameter. Typically, in ferromagnetic metals the exchange field is 
considerably higher than the internal magnetic field and it dominates the properties of the system. However, 
in some cases, where both fields can have opposite directions, an effective compensation of the conduction 
electrons polarization can occur and consequently superconductivity can be recovered at high fields H^xt — —H^x 
(Jaccarino and Peter [TD]). Bulaevskii et al. [TT] gave an excellent overview of both experimental and theoretical 
aspects of coexistence of superconductivity and ferromagnetism where both orbital and exchange effects are 
taken into account. 

The progressive development of material deposition techniques and the advent of refined lithographic 
methods have made it possible to fabricate superconductor-ferromagnet structures (S/F) at nanometer scales. 
Unlike the investigations dealing with the coexistence of superconductivity and ferromagnetism in ferromagnetic 
superconductors (for review see Flouquet and Buzdin [12|). the ferromagnetic and supercondicting subsystems in 
the artificial heterostructures can be physically separated. As a consequence, the strong exchange interaction is 
limited to a certain distance around the S /F interface whereas the weaker electromagnetic interaction can persist 
to longer distances into each subsystem. In recent reviews, Izyumov et al. |13| . Buzdin |14| and Bergeret et al. 
[T5] discussed in detail the role of proximity effects in S/F hete;ostructures dominated by "change interactionfl 
In order to unveil the effect of electromagnetic coupling it is imperative to suppress proximity effects by 
introducing an insulating buffer material between the S and F films. In an earlier report, Lyuksyutov and 

t In particular, trilayercd S/F/S structures with transparent S/F interfaces allow to realize Josephson junctions with an arbitrary 
phase difference between the superconducting electrodes, which depend on the thickness of the ferromagnetic layer (see, e.g., papers 
of Prokic et al. [16j . Ryazanov et al. [17) . Kontos et al. [18) . Buzdin and Baladie )19) . Oboznov et al. [20[ and references therein). 
The antipode F/S/F heterostructures attract a considerable attention in connection with the investigation of unusual properties of 
such layered hybrid structures governed by the mutual orientation of the vectors of the magnetization in the "top" and "bottom" 
ferromagnetic layers (see, e.g., papers of Deutscher and Meunier 1211 . Ledvij et al. |22| . Buzdin et al. I23| . Tagirov 1241 . Baladie et 
al. Gn et al. [26l|27], Pena et al. [28], Moraru et al. [29], Rusanov et al. [30], Steiner and Ziemann [31], Singh et al. [32]). 
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Figure 1. (color online) Typical examples of considered S/F hybrid systems with dominant orbital interaction. 
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Flux-coupled S/F hybrids: Dynamics of research activity 
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Figure 2. (color online) The histogram shows an increase of the number of publications dealing with the 
investigations of the S/F hybrids where the conventional low— Tc superconductors interact with magnetic 
textures mainly via stray magnetic fields: blue bars correspond to the experimental papers, while white bars 
refer to pure theoretical contributions. 

Pokrovsky [33] addressed the physical imphcations of both electromagnetic coupling and exchange interaction 
in the S/F systems deep into the superconducting state. 

In the present review we are aiming to discuss the thermodynamic and low-frequency transport phenomena 
in the S/F hybrid structures dominated by electromagnetic interactions. We focus only on the S/F hybrids 
consisting of conventional low— Tc superconductors without weak linkj[j|. The S/F heterostructures with pure 
electromagnetic coupling can be described phenomenologically using Ginzburg-Landau and London formalisms 
rather than sophisticated microscopical models. Some typical examples of such structures found in the literature 
are shown schematically in Fig. [1] As an illustration of the continuous growth of interest in S /F heterostructures 
with suppressed proximity effect we refer to Fig. [21 which shows the number of the publications during the last 
two decades. 

The review is organized as follows. Section [2] is devoted to the nucleation of the superconducting order 
parameter under inhomogeneous magnetic fields, induced by single domain walls and periodic domain structures 
in plain ferromagnetic films or by magnetic dots. A similar problem for individual symmetric microstructures 
was reviewed by Chibotaru et al. [40j. Section [3] is devoted to the static and dynamic properties of S/F systems 
at low temperatures when the superconducting OP becomes fully developed and the screening effects cannot 
be disregarded any longer. The vortex pinning properties of the S/F hybrids have been recently analyzed by 
Velez et al. [41] and the fabrication of ordered magnetic nanostructures has been earlier considered by Martin 
et al. [4^ . In the last section |4] we briefly introduce the problem of "soft" magnets in combination with 
superconducting materials, where now the superconducting currents and the magnetic stray field emanating 
from the ferromagnetic material mutually influence each other. In the conclusion, we formulate a number of 
relevant issues that, to our understanding, remain unsettled and deserve further investigations. The appendix 
summarizes the experimental and theoretical research activities on the considered S/F heterostructures, where 
we present a classification based on the choice of materials for the experimental research and on the used model 
for the theoretical treatment. 

Importantly, we would like to note already in the Introduction that the literature and references used 
by the authors in this review by no mean can be considered as a complete set. Due to dynamic and rather 
complex character of the subject and also to the limited space in this review, inevitably quite a lot of important 

X Ferromagnetic dots are shown to induce an additional phase difference in Josephson junctions, leading to a significant modification 
of the dependence of the Josephson critical current /c on the external magnetic field H^^t (so-called Fraunhofer diffraction pattern, 
see, e.g., textbook of Barone and Paterno 1341 ). which becomes sensitive to the magnetization of ferromagnetic particle (Aladyshkin 
et al. 1351 . Vdovichev et al. 136] , Fraerman et al. [371 , Held et al. 1381 . Samokhvalov '39'). 
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and interesting contributions could have been missed and, therefore, in a way, the used references reflect the 
"working list" of publications the authors of this review are dealing with. 
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2. Nucleation of superconductivity in S/F hybrids (high-temperature Umit) 

2.1. Ginzhurg- Landau description of a magnetically coupled S/F hybrid system 
Derivation of the Ginzhurg- Landau equations 

In order to describe hybrid structures, consisting of a type-II superconductor and a ferromagnet, for the case that 
no diffusion of Cooper pairs from superconductor to ferromagnet takes place, the phenomenological Ginzburg- 
Landau (GL) theory can be used. As a starting point we consider the properties of S/F hybrids for external 
magnetic fields Uext below the coercive field of the ferromagnet which is assumed to be relatively large. In this 
case the magnetization of the ferromagnet M is determined by the magnetic history only and it does neither 
depend on Uext nor on the distribution of the screening currents inside the superconductor. Such a "hard- 
magnet approximation" is frequently used for a theoretical treatment and it appears to be approximately valid 
for most part of the experimental studies presented in this section. The review of the properties of hybrid S/F 
systems consisting of superconductors and soft magnets will be presented later on in section U 

Following Landau's idea of phase transitions of the second kind, the equilibrium properties of a system close 
to the phase transition line can be obtained by minimization of the free energy functional (see, e.g, textbooks 
of Abrikosov gg, Schmidt glj, Tinkham gS]): 

II /p ^ 




V 



2 4m 



2e 

-ihV^ A^/ 

c 



.^-B.M-^},V. (1) 

where the integration should be performed over the entire spac^. Here Gso is a field- and temperature- 
independent part of the free energy, a = ao (Tcq — T), ao and (3 are positive temperature-independent constants, 
is an effective wave function of the Cooper pairs, B(r) = rot A(r) is the magnetic field and the corresponding 
vector potential, T^o is the critical temperature at _B = 0, e and m are charge and mass of carriers (e.g., 
electrons), and c is the speed of light. The term Gm, which will be explicitly introduced in the last section [4l 
accounts for the self-energy of the ferromagnet which depends on the particular distribution of magnetization. 
This term seems to be constant for hard ferromagnets with a fixed distribution of magnetization, therefore it does 
not influence the OP pattern and the superconducting current distribution in hard S/F hybrids. Introducing a 
dimensionless wave function ijj ~ ^'/^'o, normalized by the OP value — \/ oqITcq — T)/b in saturation, one 
can rewrite Eq. ([1]) in the following form 




2n 

V^p + i—Aip 
<I>o 



(2) 



expressed via the temperature-dependent coherence length = h/[Amao{Tco — T)], the London penetration 
depth — mc'^b/lSne'^ ao(Tco — T)], and the magnetic flux quantum $0 = Trhc/\e\. 

Although the Ginzburg-Landau model was proven to be consistent only at temperatures close to the 
superconducting critical temperature (Gorkov [15]), the applicability of this model seems to be much broader, 
at least from a qualitative point of view. After minimization of the free energy functional Eq. ^ with respect to 
the order parameter (OP) wave function ^ and A respectively, one can derive the two coupled Ginzburg-Landau 
equations [151 [M US] : 

-e'(v + *^A^ V'-V'+I^P^ = 0, (3) 
An 4:TT 

rot rot A = — js + 47r rot M -f — jext, (4) 
c c 

t Hereafter we used the Gauss (centimeter-gram-second) system of units, therefore all vectors B, M, and Uext have the same 
dimensionality: [_B]=Gauss (G), [M]=Oersted (Oe), [Hext]=Oersted (Oe). 
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where 



ve - 

represents the density of superconducting currents, while iext = (c/47r) rot Hea;t is the density of the currents 
corresponding to external sources, and 6 is the OP phase, ^p{r) = /(r) e*®^''^. 

Linearized GL equation 

It is quite natural to expect that at the initial stage of the formation of superconductivity [i.e. close to the 
phase transition line Tc{Hext), which separates the normal and superconducting state in the T — H^xt plane], the 
density of the superconducting condensate will be much smaller than the fully developed OP value: jV'P ^ 1- 
This allows one to neglect: (i) the nonlinear term in Eq. ([3]) and, (ii) the corrections to the vector potential 

A caused by the screening currents in Eq. (jl]), since the supercurrents js are also proportional to ji/'P- Thus, 
the nucleation of superconductivity can be analyzed in the framework of the linearized GL equation [43|, I44|, I45j : 

V + i^Aj V'^^V', (5) 
in a given magnetic field described by the vector potential distribution 

The solution of Eq. fH) consists of a set of eigenvalues (l/^^)n, corresponding to the appearance of certain OP 
pattern ipn , for every value of the applied magnetic field Hf,xt ■ The critical temperature of the superconducting 
transition Tc is determined by the lowest eigenvalue of the problem: Tc = Tco {l — {^l^^)min} ■ 

The phase boundary for plain superconducting films 

First, we would like to present the well-known solution of the linearized GL equation Eq. ([5]) corresponding to the 
OP nucleation in a plain superconducting film, infinite in the lateral direction and placed in a transverse uniform 
magnetic field 'Rext = H^xt Zq [131 Sll HSj ■ Taking the gauge Ay = xH^xt : one can see that Eq. ^ depends 
explicitly on the x— coordinate only, therefore its general solution can be written in the form: ip — f(x) e''=f+"?^, 
where the wave vectors k and q should adjust themselves to provide the maximization of the Tc value. Using 
this representation in Eq. ([5]), it is easy to see that the spectrum of eigenvalues (1/C^)n is similar to the energy 
spectrum of the harmonic oscillator but shifted: (1/C^)n = 27r(2n + l)Hext/^o + q^ and the (1/C^) minimum 
(the maximum of Tc) corresponds to n = and g = for any Hext value, (l/f^)min = '^irHcxt/^o- The critical 
temperature of the superconducting transitiorQ as a function of a uniform transverse magnetic field is given by 
Tc = Tco [1 - 2<2 Hcxt/'i>o] , or 

where H^^ = $o/(27r^o) is the upper critical field at T = 0. The inversely proportional dependence of the 
shift of the critical temperature 1 — Tc/Tco on the square of the OP width ifj can be interpreted in terms 
of the quantum-size effect for Cooper pairs in a uniform magnetic field. It should be mentioned that the 
effect of sample's topology on the eigenenergy spectrum {l/£,'^)n becomes extremely important for mesoscopic 
superconducting systems, whose lateral dimensions are comparable with the coherence length ^. Indeed, 
this additional confinement of the OP wave function significantly modifies the OP nucleation in mesoscopic 
superconductors and the corresponding phase boundaries Tc{Hext) differ considerably from that typical for 
bulk samples and films infinite in the lateral directions (see Chibotaru et al. [40] , Moshchalkov et al. 07[ [48] , 
Berger and Rubinstein 49j). 

t It is well known that superconductivity nucleates in the form of a Gaussian-Uke OP wave function ip{x,y) = e~^^~^o' /^^jf e''°^, 
locahzed in the lateral direction at distances of the order of the so-called magnetic length fi^ = $o/(27r|i?ea;t |) and uniform over 
the film thickness. The oscillatory factor e'*^ describes the displacement of the OP maximum positioned at xq = k^Q/{2-KH^xt) 
without a change of the (l/5^)min value. It is interesting to note that the confinement of the OP wave function is determined by 
the magnetic length lu, i.e. the OP width is a function of the external field Hext- On the other hand, the temperature-dependent 
coherence length ^ is a natural length scale describing the spatial OP variations. The equality = defines the same phase 
boundary in the T — Hext plane as that given by Eq. (71 . 
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2.2. Magnetic confinement of the OP wave function in an inhomogeneous magnetic field: general 
considerations 

The main focus in this section is to describe the nucleation of the superconducting order parameter in a static 
non-uniform magnetic field He^t + b(r) based on a simple approaclQ- This method makes it possible to see 
directly a correspondence between the position of the maximum of the localized wave function -0 and the 
critical temperature Tc in the presence of spatially-modulated magnetic field b(r), generated by ferromagnet. 
For simplicity, we assume that the thin superconducting film is infinite in the (x, j/)— plane, i.e. perpendicular 
to the direction of the external field He^t = Hext'^o- This allows us to neglect the possible appearance of 
superconductivity in the sample perimeter (surface superconductivity [1211311133]) and focus only on the effect 
arising from the nonuniform magnetic field. 

Importance of out- of -plane component of the field 

It should be emphasized that the formation (or destruction) of superconductivity in thin superconducting films 
is sensitive to the spatial variation of the out-of-plane component of the total magnetic field. Indeed, the upper 
critical fields and ijjj for the out-of-plane and in-plane orientation for a uniform applied magnetic field can 
be estimated as follows gH |44l |45] : 

r± ^ rrW ^ 

where Ds is the thickness of the superconducting sample. For rather thin superconducting films and/or close 
to the superconducting critical temperature Ds ^ i — T /Tcq)^^/^ , therefore <C h\2- other 

words, superconductivity will generally be destroyed by the out-of-plane component of the magnetic field 
rather than by the in-plane component, and thus, to a large extent, the spatial distribution of the out-of- 
plane component determines the OP nucleation in thin-film structures. Since a uniform magnetic field is known 
to suppress the critical temperature, one can expect that the highest Tc value should correspond to the OP 
wave function localized near regions with the lowest values of the perpendicular magnetic field |i?2(r)| provided 
that = Hext + ^z(r) varies slowly in space. 

If the field Hgxt exceeds the amplitude of the internal field modulation (i.e. Hext < —uiaxbz and 
Hext > ^miii6^), the total magnetic field is non-zero in the whole sample volume, and the favorable positions 
for the OP nucleation are at the locations of minima of |i?z(r)| = \Hext + ^z(r)|. If the characteristic width 
of the OP wave function, which will be defined later, is much less than the typical length scale £b of the magnetic 
field variation, then locally the magnetic field can be considered as uniform at distances of the order of and 
it approximately equals to iiim\Hext + &z(r)|. Then, using the standard expression for the upper critical field 
Eq. ([7|) and substituting the effective magnetic field instead of the applied field, one can obtain the following 
estimate for the phase boundary 

1 - — ~ mm , li/eitl > max lO^I. (9j 

According to this expression, the dependence Tc{Hext) is still linear asymptotically even in the presence of a 
nonuniform magnetic field. However, the critical field will be shifted upwards (for H^xt > 0) and downward (for 
Hext < 0) by an amount close to the amplitude of the field modulation. Generally speaking, such a "magnetic 
bias" can be asymmetric with respect to the Hext = provided that max5z(a;) ^ \mm.hz{x)\. 

For relatively low Hext values, when the absolute value of the external field is less than the amplitude of 
the field modulation (— maxfo^ < Hext < — rninfe^), the z— component of the total magnetic field Hext + bz{x) 
becomes zero locally somewhere inside the superconducting film. As a result, superconductivity is expected to 
appear first near the positions where Hext + ^z(ro) = [see panel (a) in Fig. [5]. It is natural to expect that 
the details of the OP nucleation depend strongly on the exact topology of the stray field as well as the field 
gradient near the lines of zero field. As an example we will analyze the formation of superconductivity in a thin 
superconducting film placed in an inhomogeneous magnetic field modulated along certain direction. 

OP nucleation in a magnetic field modulated in one direction 

t We introduce the following notations: b = rota characterizes the non-uniform component of the magnetic field only, while 
B = Hext + b = rot A is the total magnetic field distribution, the external field He^^t is assumed to be uniform. 
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Figure 3. (color online) (a) Schematic representation of the OP wave function tl){x) localized near the point 
xo, where the z— component of the total magnetic field Bz = H^^t + bz vanishes. Provided that the OP width is 
much smaller than the typical length scales of the magnetic field <C £(,), the actual field distribution Bz(x) 
can be approximated by a linear dependence B^ix) ~ (dbz/dx)^^ (x — xq). 
(b). Energy spectrum eo vs. Q of the model problem Eq. Illll l. after Aladyshkin et al. |51| . 

Following Aladyshkin et al. [SOj, we estimated the dependence of Tc{Hext) for a thin superconducting film 
in the presence of a non-uniform magnetic field modulated along the a;— direction, where (x,y,z) is the 
Cartesian reference system. Let the external field be oriented perpendicular to the plane of the superconducting 
film, Hext = H^xt'^Oi while the z— component of the total magnetic field vanishes at the point xq, i.e. 
Hext + bz{xo) — 0. The vector potential, corresponding to the field distribution Bz{x) — Hext + bz{x), can 
be chosen in the form Ay(x) = xH^xt + Oy(a^), where bz = duy/dx. Since there is no explicit dependence on 
the J/— coordinate in the linearized GL equation Eq. ([5]), the solution uniform over the sample thickness can be 
generally found as ip{x, y) = fk{x) e**^^ and the absolute value of the OP satisfies the following equation: 

Now the parameter k cannot be excluded by a shift of the origin of the reference system, therefore one should 
determine the particular k value in order to minimize (l/^o)(l — Tc/Tco) and thus, to maximize the Tc value. 

For an unidirectional modulation of the field, the curves of zero field, where we expect the preferable 
OP nucleation, are straight lines parallel to the y— axis, and their positions depend on the external field, 
xo — xo{Hext)- Expanding the vector potential inside the superconducting film in a power series around the 
point Xo, one can get 

Ay{x) ~ xoHext + ay{xo) + ^ b'^{xo)(x ~ Xq)"^ + ... 

This local approximation is valid as long as |6"(a;o)^^/6^(xo)| ^ 1. Introducing a new coordinate r = {x—xo)/£jp 
and the following auxiliary parameters £^ and Qk 

we can reduce Eq. (jlOp to the bi-quadratic dimcnsionless equation 

Thus, the problem of the calculation of the highest Tc value in the presence of an arbitrary slowly-varying 
magnetic field as a function of both the external field and the parameters of the "internal" fields is reduced 
to the determination of the lowest eigenvalue Eq — So{Q) of the model equation Eq. pT|) . As was shown in 
Ref. [51], the function eo{Q) is characterized by the following asymptotical behavior: eo{Q) ~ -I- \J—Q for 
Q <C — 1 and eo(Q) — for Q » 1 and it has the minimum value e^nin = 0.904 [the panel (b) in Fig. [3]. 
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Extracting T^ from £q = (1 — Tc/Tcq) ■ ^^/Coj the approximate expression of the phase boundary takes the 
foUowing form: 

C2 / I?,/ ('^„M \ 2/3 



I- mineo(Qfe) - f^^^Nr^) . -max b, < H,,t < -min b,. (12) 



If there are several points xo,i where the external field compensates the field generated by the ferromagnetic 
structure, then the right-hand part of Eq. (fT2)) should be minimized with respect to xo^i. The application of 
Eg. p^ for the model cases bz = 4Ms arctan(Z?//a;) (single domain wall) and bz = Bq sm{2Trx / w) (periodic 
domain structure) were given in Ref. [50j . 

OP nucleation in axially- symmetrical magnetic field 

Similar to the discussion above, one can expect that in the presence of an axially-symmetrical magnetic field 
superconductivity will nucleate in the form of ring-shaped channels of radius r = rg, where H^^xt + ^z(^o) ~ 0- 
The independence of the linearized GL equation Eq. ^ on the angular (yS— coordinate results in a conservation 
of the angular momentum (vorticity) L of the superconducting wave function. Thus, a nonuniform magnetic 
field makes it possible to have an appearance of giant (multiquanta) vortex states, which are energetically 
unfavorable in plain (non-perforated) large-area superconducting films, but have been observed in mesoscopic 
superconductors (Moshchalkov et al. [47l|48], Berger and Rubinstein [49]) and nanostructured films with antidot 
lattices (Baert et al. [52j, Moshchalkov et al. [53]). Expanding the vector potential in the vicinity of tq and 
repeating similar transformations as above, one can get the following approximate expressioijil for the phase 
transition line (Aladyshkin et al. [STl I54j ) 



1 - — ^ ~ mineo(QL) - j^- ^ -max bz < H^xt < -min bz 



where the parameters £^ — y^o/{n \dbz/dr\^^^) and Ql = — [27rroyl;^(ro)/$o — L] £^/ro depend on the external 
field. Thus, this model predicts field-induced transitions between giant vortex states with different vorticities. 
Since the vorticity L is a discrete parameter, the changes of the favorable L value while sweeping the external 
field leads to abrupt changes in dTc/dHext- Similar periodic oscillations of Tc were originally observed by Little 
and Parks |551 I56j for a superconducting cylinder in a parallel magnetic field and later for any mesoscopic 
superconductor in a perpendicular magnetic field (for a review see Chibotaru et al. [40j ) . 

Effect of nonuniform magnetic field on two-dimensional electron gas 

It is interesting to note that there is a formal similarity between the linearized GL equation ^ and the stationary 
Schrodinger equation (see, e.g., Ref. fS7j) for a charged free spinless particle in a magnetic field 



2 / ■ \ 2 

V-^A ^ = Ei^, (13) 



h 

2m V he 



where "0 is the single particle wave function, e is the charge and m is the mass of this particle. Based on this 
analogy, one can map the results, obtained for a normal electronic gas in the ground state, on the properties of 
the superconducting condensate near the "superconductor-normal metal" transition on the T — H^xt diagram. 

In particular, the effect of a unidirectional magnetic field modulation on the energy spectrum of a two- 
dimensional electronic gas (2DEG) was analyzed by Miiller [35], Xue and Xiao ^55], Peeters and Vasilopoulos 
[50] . Peeters and Matuhs [SI], Wu and UUoa [S5], Matulis et al. [S3], Ibrahim and Peeters [M], Peeters et al. 
[55] . Gumbs and Zhang [55], Reijniers and Peeters [S71|5H]j Nogaret et al. [MllTD]. Remarkably, in a magnetic 
field varying linearly along certain direction, quasi-classical electronic trajectories propagating perpendicularly 
to the field gradient V-B^ are confined to a narrow one-dimensional channel localized around the region where 
Bz = [58j . A more detailed numerical treatment revealed a lifting of the well-known degeneracy of the Landau 
states on the centers of the Larmor's orbit, inherent to electrons in a uniform magnetic field. Periodic magnetic 
field patterns with zero and nonzero average value was shown to transform the standard Landau spectrum 
E = fiujc{n + 1/2) into a periodic E{ky) dependence, describing the broadening of the discrete Landau levels 

t Note that the case of a point magnetic dipole approximation fails for L = since the maximum of the OP wave function is 
located at r = and this theory cannot correctly describe neither the OP nucleation nor the phase boundary Tc(Hext) at negative 
Hext values close to the compensation field Bq = — ma'Kbz{r). 
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into minibands as in the case of one-dimensional potential (here lOc = \e\H(.,j.t/ {mc) is the cyelotron frequency). 
An oscillatory change of the width of energy bands as H is swept was shown to give rise to oscillations in the 
magnetoresistance of 2DEG at low H^xt values, which reflect the commensurability between the diameter of 
cyclotron orbit at the Fermi level and the period of the magnetic field modulation [55] • The mentioned oscillatory 
magnetoresistance due to commensurability effects was later on corroborated experimentally by Carmona et 
al. |71| . The influence of two-dimensional magnetic modulations on the single-particle energy spectrum was 
theoretically considered by Hofstadter J72j. A modification of the scattering of two-dimension electrons due to 
the presence of either ferromagnetic dots or superconducting vortices was shown to lead to a nontrivial change 
of the conductivity in various hybrid systems: 2DEG/superconductor (Geim et al. [73 , Brey and Fertig [74] . 
Nielsen and Hedegard Tl], Reijniers et al. [W) and 2DEG/ferromagnet (Khveshchenko and Meshkov [77], Ye 
et al. ,78j, Solimany and Kramer [79j . Ibrahim et al. 80J, Sim et al. 81j, Dubonos et al. |82| . Reijniers et 
al. [551 184j ). The effect of inhomogeneous magnetic field on the weak-localization corrections to the classical 
conductivity of disordered 2DEG was considered by Rammer and Shelankov [85] , Bending [86] , Bending et al. 
[571IHH1IMI, Mancoff et al. [SD], Shelankov [H], Wang 



2.3. Planar S/F hybrids with ferromagnetic bubble domains: theory 

The aforementioned Ginzburg-Landau formalism can be applied in the case of a non-uniform magnetic field 
generated by the domain structure of a plain magnetic film. The problem of the OP nucleation in planar 
S/F hybrid structures was theoretically analyzed for hard ferromagnets characterized by an out-of-plane 
magnetization M = (x) Zq by Aladyshkin et al. [50 , Buzdin and Melnikov [Q^ , Samokhin and Shirokoff ^4] , 
Aladyshkin et al. [95], Gillijns et al. [96 . It is also worth mentioning the pioneer paper of Pannetier et al. [97] . 
where the OP nucleation in a periodic sinusoidal magnetic field generated by a meander-like lithographically- 
prepared metallic coil was considered. 

The distribution of the vector potential a(r) can be obtained, either by integration of the last term in the 
r.h.s. of Eq. ([5]) or by a direct consideration of the magnetostatic problem. Provided that the width of the domain 
walls S is much smaller than other relevant length scales, the field distributions can be calculated analytically for 
some simple configurations (Aladyshkin et al. [55], Sonin [551115]). Choosing the gauge Ay = xHext + ay{x, z), 
one can easily see that the linearized GL equation Eq. ([5]) does not depend on the y— coordinate, hence we 
can generally find the solution in the form i^i^: Vi ^) — fk{x, z) exp(zfcy), where the function /fc(a;, z) should be 
determined from the following 2D equation: 



, , 27r 
—ay{x,z) + —xHext 
5>0 "Po 



/fe = ^/fc, (14) 



If the superconducting film has insulating interfaces at the top and bottom surfaces and A • n = at the surface 
dVs of the superconductor, one should apply the standard boundary conditions: dfk/dn = (here n is the 
normal vector). 

The spatial distribution of the magnetic field, induced by the periodic ID domain structure, strongly 
depends on the relationship between the width of the magnetic domains w and the thickness of the ferromagnetic 
film Df as well as on the distance h between superconductor and ferromagnet (Fig. [J). If the superconducting 
film thickness Dg is much smaller than the typical length scales of the nonuniform magnetic field {w and Df), 
in a first approximation one can neglect the OP variations in the z— direction and omit the term d^fk/dz^ in 
Eq. (fT4|) . As a result, the OP nucleation in a thin superconducting film is determined by the spatial profile of 
the perpendicular magnetic field only, while the effect of the parallel field can be ignored. 

Criterium for the development of domain-wall superconductivity 

In this section we will discuss the possibility of localizing the OP wave function near a domain wall at zero 
external field. Obviously, the regime of domain-wall superconductivity (DWS) can be achieved only if the 
typical 5* value inside the magnetic domains is rather large in order to provide an exponential decay of the 
OP, described by the effective magnetic length il = ^/Wo/2^T\b*\ , within a half-width of the domain: £1 < w/2. 
For thick ferromagnetic films {w/Df <C 1) the magnetic field inside domains is almost uniform and it can be 
estimated as 5* ~ 2ttMs, giving us the rough criterion of the realization of the DWS regime and the critical 
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Figure 4. (Color online) Transverse 2— component of the magnetic field, induced by one-dimensional periodic 
distribution of magnetization with the amplitude Ms and the period w, calculated at the distance h <^ Df 
above the ferromagnetic film of a thickness D j , adapted from Aladyshkin and Moshchalkov |95| . 

temperature TjP^ at H^xt = 0: tt^A/^wV^o > 1 and (T^o - Tj:°^)/Tco 2'kMs/H^^\ Of course, 2t:M,/H^°^ 
should be less than unity otherwise superconductivity will be totally suppressed. 

By applying an external field on the order of the compensation field, Hext — SttM^, one can get local 
compensation of the field above the domains with opposite polarity and a doubling of the field above the 
domains of the same polarity. Since superconductivity is expected to form at regions with zero field (which 
are w wide), the maximal critical temperature can be estimated as follows: (Tco — T™°-^)/Tco — Co/^^ 
consequence of the quantum-size effect for Cooper pairs in nonuniform magnetic field). Therefore, T™"^ will 




H,,/(4M^) w/D, 

Figure 5. (Color online) (a) Examples of the phase transition lines Tc{Hext) for a planar S/F structure 
containing a periodic ID domain structure (Ms is the saturated magnetization, Df is the ferromagnetic film 
thickness, w is the period of the domain structure, the thickness of the superconducting film Ds ^ {Df,w) 
and the separation between superconducting and ferromagnetic films h/ Dj <^ 1), adapted from Aladyshkin and 
Moshchalkov [95] . Black dashed line correspond to the Tc(Hext) dependence in the absence of the non- uniform 
field. 

(b) Different regimes of localized superconductivity in the presence of a ID domain structure in the Ms — w 
plane, obtained numerically for Ds <C {Df,w) and h/Df <^ 1, adapted from Aladyshkin and Moshchalkov 1951 . 
In regions II and III the phase boundary Tc{Hext) exhibits reentrant superconductivity. The slope dTc/d\Hext\ 
at Hext = can be positive (III), zero (II) or negative (II, near the separating line I— II). Region I corresponds 
to the monotonic Tc{Hext) dependence. 
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exceed Tc , pointing out to the non-monotonous Tc{Hext) dependence for the same Alg and w parameters, 
which are necessary to have the DWS regime at H^xt = 0. The typical phase boundary Tc{Hext), corresponding 
to the DWS regime at Hext = and shown in panel (a) in Fig.O is characterized by the presence of a pronounced 
reentrant behavior and the parabolic dependence of Tc on Hext at low fields (curve labelled w/Df = 8). This 
type of phase boundary was predicted by Pannetier et al. j97j for a superconducing film in a field of parallel 
metallic wires carrying a dc current, and by Buzdin and Melnikov [93] and Aladyshkin et al. ^501 195j for planar 
S/F hybrids. 

Localized superconductivity in S/F hybrids for w/Df <^ 1 and w/Df ^ 1 

For S/F hybrids with smaller periods of the field modulation {tt^MsW^ /^o ^ 1) the OP distribution cannot 
follow the rapid field variations and, as a consequence, at H^xt — there is a broad OP wave function, spreading 
over several domains and resulting in an effective averaging of the nonuniform magnetic field. In this case the 
critical temperature was shown to decrease monotonically with increasing |-ffea:t| [curve labelled w/Df = 2.5 in 
Fig. Elja)], similar to the case of superconducting films in a uniform magnetic field. By applying an external 
field, one can shrink the width of the OP wave function and localize it within one half-period above the 
domains with opposite magnetization. The interplay between both, the external field and the periodic magnetic 
field, which determines the resulting OP width, leads to a sign change of the second derivative d^Tc/dH^^^. 
At high Hcxt values the width of the OP wave function, positioned at the center of the magnetic domain, 
is determined by the local field Bioc — {H^xtl — SttM^, therefore we come to a biased linear dependence 
1 — Tc/Tco ~ jl-ffeitl — 27r Ms I /iJ|^2'' • These qualitative arguments were supported by numerical solutions of 
the linearized GL equation [95] . 

The case w/Df ^ 1 should be treated separately since the z— component of the field inside the magnetic 
domains is very inhomogeneous: the absolute value |62(a;)| reaches a minimum 6* = SirMgD f /w at the domain 
center, while the maximal value is still equal to 2'kAIs at the domain walls. It was shown that the |?>2(a^)| 
minima are favorable for the OP nucleation at Hext = 0. In this regime, the OP localization in the center of the 
domains at Hgxt = is possible as long as 2-K^MswDf /^q > 1. At the same time the nucleation near domain 
walls is suppressed by the mentioned field enhancement near the domain walls. The sudden displacement of the 
localized OP wave function between the centers of the domains of positive and negative magnetization, when 
inverting the Hext polarity, results in a new type of phase boundary Te{Hext) with a singularity at Hext = 
|95j . It is important to note that, for w/Df ^ 1 and Hext = the critical temperature increases linearly with 
almost the same slope dTc/ d\Hext\ — Teo/H^ as the Tc value decreases in an applied uniform magnetic field 
[curve labelled w/Df = 14 in Fig. [5] (a)]. 

2. 4- Planar S/F hybrids with ferromagnetic bubble domains: experiments 
OP nucleation in perpendicular magnetic field 

To the best of our knowledge, the first observation of reentrant superconductivitjQ in planar S /F hybrids was 
reported by Yang et al. |101j who measured the electrical resistance of a superconducting Nb film grown on top 
of a ferromagnetic BaFei20i9 substrate characterized by an out-of-plane magnetization. Later on, the same 
system Nb/BaFei20i9 was examined by Yang et al. in Ref. [102j . From the parameters typical for the domain 
structure in BaFei20i9 single crystals and Nb films {Mg ~ 10^ Oe, w ~ 2 /im, Df ~ 90 /xm, iJ^j"* — 30 kOe), the 
following estimates can be obtained: w/Df ~ 0.02, tt^M^wV^o > 10^ and 2T:MjH^f ~ 0.02. Therefore, such 
a ferromagnet is suitable for the realization of the DWS regime at Hext = 0. The appearance of these localized 
superconducting paths guided by domain walls was shown to result in a broadening of the superconducting 
resistive transition at low magnetic fields. As the field Hext is ramped up, the superconducting areas shift away 
from the domain walls towards the wider regions above the domains with an opposite polarity (so-called reversed- 
domain superconductivity, RDS) where the absolute value of the total magnetic field is minimal because of the 

t The experimental observation of tiie influence of a periodic magnetic field, generated by an array of parallel wires with current 
/ flowing alternatively in opposite directions, on the properties of an Al superconducting bridge was reported by Pannetier et al. 
|97j . Since the max|6z| oc /, reentrant superconductivity can be realized for rather high / values, as it was shown experimentally. It 
should be noted that already in sixties Artley et al. llOOl experimentally studied the effect of the domain walls in a thin permalloy 
fllm on the superconducting transition of a thin indium fllm. 
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Figure 6. (color online) Preparation of the magnetic state in a ferromagnetic Co/Pt film with a desirable 
remanent magnetization Mrem- 

(a) Magnetization loops M{Hext) at 300 K (triangles) and 5 K (squares), the magnetic field axis is normalized 
by the corresponding coercive fields H^, adapted from Gillijns et al. |96h 

(b) Remanent magnetization Mrem, measured at 5 K and Hext = after saturation in positive fields (up to 10"* 
Oe) and subsequent application of a returning field Hret [this procedure is shown schematically in panel (c)], 
adapted from Gillijns et al. I9bl . 



compensation effect. As a consequence, the superconducting critical temperature Tc increased with increasing 
l-ffeitl up to 5 kOe. Once the external field exceeds the saturation field Hg of the ferromagnet {Hg ~ 5.5 kOe 
at low temperatures), the domain structure in the ferromagnet disappears and the phase boundary abruptly 
returns back to the standard linear dependence (1 — Tc/Tco) — \H),,j:t\/ H^^\ Since the width and the shape 
of the magnetic domains continuously depend on the external field, the theory developed in section 12.21 is not 
directly applicable for the description of the experiment, although it qualitatively explains the main features of 
the OP nucleation in such S/F systems. 

Substituting Nb by a superconductor with a smaller iJj^2'' value (e.g., Pb with -ff^^'* — 1-8 kOe) allows one 

to study the effect of the superconducting coherence length — \J '^o/'^t^H^ on the localization of the OP. It 
was shown by Yang et al. [103] that the increase of the Ms/h'^ ratio suppresses the critical temperature of the 
formation of domain- wall superconductivity at zero external fielcjj, therefore superconductivity in Pb/BaFei20i9 
hybrids appeared only near the compensation fields above the reversed domains. 

Direct visualization of localized superconductivity in Nb/PbFei20i9 structures was performed by Fritzsche 
et al. |104j . The basic idea of this technique is the following: if the sample temperature becomes close to a 
local critical temperature at a certain position then a laser pulse, focused on that point, can induce the 

local destruction of superconductivity due to heating. The observed increase of the global resistance R of the 
superconducting bridge can be associated with the derivative dR{x,y)/dT. By varying the temperature and 
scanning the laser beam over the Nb bridge under investigation, it is possible to image the areas with different 
critical temperatures. For example, it allows one to attribute the formation of well-defined regions with rather 
high local critical temperatures above magnetic domains at the compensation field with the appearance of 
reversed-domain superconductivity. 

The effect of the amplitude of the field modulation on the OP nucleation was considered by Gillijns et al. 
[10511106] on thin-film trilayered hybrid F/S/F structures. In contrast to the BaFei20i9 single crystal discussed 
above, the multilayered Co/Pd films are characterized by a high residual out-of-plane magnetization, Mg ~ 10^ 
Oe, almost independent on the external field at |i?ea:t| < Hcoer, where Hcoer — 10"^ Oe is the typical coercive 
field at low temperatures. The use of two ferromagnetic films with slightly different coercive fields allowed 
them to prepare different magnetic configurations and thus to control the amplitude of the nonuniform field 

I In section [2. 31 we argued that the critical temperature, tJ"', at zero external field is proportional to f — 2iTMa / . 
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Figure 7. (color online) (a)-(d) MFM images obtained at T = 300 K for Hret values equal to -1.75 kOe (a), 
-2.00 kOc (b), -2.50 kOe (c), -3.00 kOe (d), the coercive field iT^OO^ = 1.91 kOe, adapted from Gillijns et al. 
1961 . The dark (bright) color represents domains with positive (negative) magnetization. 

(e) A set of experimental phase boundaries Tc{Hext) obtained for the same bilayered S/F sample (a 
superconducting Al film on top of a Co/Pt multilayer) in various magnetic states measured after the procedure 
of an incomplete demagnetization: Hext = Hext = 10 kOe =^ Hext = Hret => Hext = for various 
returning fields Hret indicated on the diagram, the coercive field H^^ = 3.97 kOe, adapted from Gillijns et al. 



inside the superconductor due to the superposition of the partial stray fields via an appropriate demagnetizing 
procedure. The effective doubling of the amplitude of the internal field for a configuration with two demagnetized 
ferromagnetic films (containing bubble domains) leads to a broadening of the temperature interval, where the 
Tc{Hext) hne demonstrates the non-monotonous behavior. In other words, the critical temperature of domain- 
wall superconductivity at H^xt — expectedly decreases as the effective magnetization increases. In addition, 
the enhancement of the internal field results in a shift of the T^ maxima to higher i/e^t values, what corresponds 
to reversed-domain superconductivity. 

The spatial extension where the field compensation takes place, is a crucial parameter defining the nucleation 
of superconductivity: an OP trapped in a broader region results in a higher Tc value and vice versa (Gillijns et 
al. [96], Aladyshkin et al. [107]). The magnetic state of the ferromagnet can be reversibly changed after the 
following procedure of an incomplete demagnetization : Hext = Hext = Hg => Hext = Hret =^ Hext = 0, 
where Hg is the saturation field (see Fig. [6]). As a result, one can obtain any desirable remanent magnetization 
—Ms < Mrern < Mg (as wcll as any average width of the magnetic domains) by varying the Hret value. At 
Hret < the formation of the negative domains decreases the average width of the positive domains and 
it causes a drastic lowering of the height of the Tc peak, positioned at negative fields and attributed to the 
appearance of superconductivity above large positive domains [curves Hret = —3.93 kOe and Hret = —4.15 kOe 
in Fig. [7|. This observation is a direct consequence of the increase of the ground energy of the "particle- in-a- 
box" with decreasing width of the box. When Mrem is close to zero, thus indicating the presence of an equal 
distribution of positive and negative domains, a nearly symmetric phase boundary with two maxima of the same 
amplitude is recovered [curve Hret = —4.55 kOe in Fig. [7]. For higher Hret values, the first peak, located at 
negative fields, disappears, whereas the peak at positive fields shifts up in temperature and is displaced to lower 
magnetic field values [curves Hret = —4.61 kOe and Hret = —5.00 kOe in Fig. [7|. This second peak eventually 
evolves into a linear phase boundary when the ferromagnetic film is fully magnetized in the negative direction. 

OP nucleation in parallel magnetic field 

We would like to mention a few related papers devoted to the nucleation of superconductivity in various planar 
S/F structures, where superconducting and ferromagnetic layers were not electrical insulated and thus an effect 
of exchange interaction can not be excluded. 

An appearance of domain walls in permalloy film leading to dips in a field dependence of the resistivity 
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of a superconducting Nb film was observed by Rusanov et al. |108| . The position of these resistivity minima 
were found to be dependent on the sweep direction of the in-plane oriented external field. An opposite effect 
(maxima of resistivity at temperatures below Tco) was observed by Bell et al. [109] for thin- film amorphous S/F 
structures consisting of superconducting MoGe films and ferromagnetic GdNi layers. It was interpreted as flow 
of weakly pinned vortices induced by the stray field of magnetic domains in the ferromagnetic layers. Zhu et al. 
[llOj demonstrated that the domain structure in multilayered CoPt films can be modified by applying in-plane 
external field during the deposition process. This deposition field does not change the overall perpendicular 
magnetic anisotropy of the Co/Pt films but it induces a weak an in-plane magnetic anisotropy and eventually 
alters the domain patterns. Indeed, after demagnetizing with an in-plane ac magnetic field oriented along the 
deposition field direction, one can prepare the domain structure in a form of largely parallel stripe domains. In 
contrast to that, the same multilayered structure fabricated at zero field or demagnetized with out-of-plane ac 
field exhibits a nearly random labyrinth- type domain pattern. Sweeping the external field ijjjjt , one can control 
the arrangement of domain walls, and drive the S/F hybrid from normal to superconducting state for the same 
temperature and magnetic field. 

2.5. S/F hybrids with 2D periodic magnetic field: theory and experiments 

In this section we continue the discussion concerning the properties of large- area S/F hybrids in which the 
nonuniform magnetic field is created by regular arrays of ferromagnetic dots. The fabrication of such magnetic 
structures makes it possible to achieve full control of the spatial characteristics of the nonuniform magnetic 
field (both the topology and the period), which can be eventually designed practically at will. One can expect 
that due to the field-compensation effect the inhomogeneous magnetic field modulated in both directions will 
affect the OP nucleation in the same way as for the planar S/F structures. However, the 2D periodicity of the 
magnetic field naturally leads to the appearance of well-defined commensurability effects for such hybrid systems. 
In other words, a resonant change in the thermodynamical and transport properties of the superconducting films 
appears as a function of the external magnetic field H^xt, similar to that observed for superconductors with 
periodic spatial modulation of their properties (e.g., perforated superconducting films [52[|53j ). These matching 
phenomena take place at particular Hext,n values that can be used as indicators allowing us to find a relationship 
between the most probable microscopic arrangement of the vortices in the periodic potential and the global 
characteristics of the considered S/F hybrids measured in the experiments. 

Commensurate solutions of the GL equations. 

The periodic solutions of the GL equations in the presence of a nonuniform 2D periodic magnetic field can 
be constructed by considering one or more magnetic unit cells (of total area S) and applying the following 
boundary conditions: 



where b^, k = {x, i/} are the lattice vectors, and rjk is the gauge potential (Doria et al. [lllj ). The gauge 
transformation Eq. ([T5)) is possible provided that the flux induced by the external magnetic fleld H^xtS through 
the chosen area S is equal to an integer number of flux quanta, n$Oi which gives us the matching fields 



The formation of different vortex patterns in S/F hybrids, containing square arrays of the magnetic dots 
with perpendicular magnetization, at H^xt — were studied by Priour and Fertig [1121 1113j and Milosevic 
and Peeters |1141 1115j . Since the total fiux through an underlying superconducting film, infinite in the lateral 
direction, equals to zero, vortices cannot nucleate without corresponding antivortices, keeping the total vorticity 
zero. It was shown that the number of the vortex-antivortex pairs depends on the dipolar moment of the 
magnetic dot. The equilibrium vortex phase, corresponding to the minimum of the free energy functional, can 
exhibit a lower symmetry than the symmetry of the nonuniform magnetic field. In order to get such vortex 
states of reduced symmetry, the GL equations with periodic boundary conditions should be considered in a 
large supercell (2 x 2, 4 x 4 etc.). In most cases vortices are confined to the dot regions, while the antivortices, 
depending on the dot's magnetization, can form a rich variety of regular lattice states with broken orientational 
and mirror symmetries (Fig. [8]). The creation of vortex-antivortex pairs in case of ferromagnetic dots with 



A(r + bfc) = A(r)+V?7fe(r), 

^'(r + bfc) = «'(r) exp (27ri77fc(r)/$o) , 



(15) 




in^o/S". 
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(b) N„=2.5, N3„=2.5 



(c) N„=3,N3„=3 



Figure 8. (color online) Contour-plots of the Cooper pair density of stable vortex phases in superconducting 
thin films in the presence of a magnetic dot array with dipolar moments oriented perpendicular to the film plane, 
adapted from Priour and Fertig |112l . The images in (a), (b), and (c) correspond to magnetic dot strengths of 
3.61, 4.33, and 4.95 fundamental flux units, respectively (the dot's dipole moment is specified by the positive 
fiux passing through each unit cell). The Cooper pair density is highest in red and lowest in blue regions. 
Antivortices appear as small blue spots between the magnetic dots depicted by black squares. The number of 
vortices and antivortices Nav per unit cell are indicated on the plots. The calculations were performed for 
a large 4x4 supercell with periodic boundary conditions and for the following parameters: k = A/^ = 4, the 
period of the magnetic dot array is 6.25 ^, the lateral size of the magnetic dot is 2^ X 2^. 



in-planc magnetization was considered by Milosevic and Peeters |116) . As expected, vortex-antivortex pairs 
appear under the poles of each magnet according to their specific inhomogeneous magnetic field keeping the 
total flux through the superconducting film equals to zero. 

The influence of the external field on the formation of symmetrical and asymmetrical commensurate vortex 
configurations in thin superconducting films in the presence of a square array of out-of-plane magnetized dots 
was investigated by Milosevic and Peeters [114[ I115[ 11171 I118j . The simulations were carried out only for 
some discrete values Hext.n of the external field, corresponding to the magnetic flux quantization per magnetic 
supercell of area S. Vortices were shown to be attracted by the magnetic dots in the parallel case (at Hext,n > 
for Mz > 0) and repelled in the antiparallel case (at Hext,n < for Mz > 0). In the parallel case the vortex 
configurations for the integer matching fields are similar to that for the vortex pinning by regular arrays of 
antidots. 

The temperature dependence of the magnetization threshold for the creation of vortex-antivortex pairs was 
considered in Ref . [118] . It was noted that the system will not necessarily relax to the ground state, if there are 
metastable states, corresponding to local minima in the free energy. As long as the given vortex state is still 
stable and it is separated from other stable vortex configurations by a finite energy barrier (analogous to the 
Bean-Livingston barrier for the vortex entry into superconducting samples), then the vorticity remains the same 
even when changing temperature. However, an increase in temperature resulting in a decrease of the height 
of the energy barriers and strengthening of the thermal fluctuations, can eventually cause a phase transition 
between the vortex states with different number of vortices. The modification of the ground state (at H^xt — 0) 
by the creation of extra vortex-antivortex pairs, when changing temperature and/or increasing the Mg value 
manifests itself as cusps in the phase boundary separating superconductor from the normal metal phase in the 
Ms — T diagram, similar to the Little-Parks oscillations in the Tc{Hext) dependence [551 156j . 

Oscillatory nature of the phase transition line (in-plane dot's magnetization) 

The influence of two-dimensional square arrays of micron-sized, in-planc magnetized particles (SmCo, GdCo, 
FeNi) on the electrical resistance of a superconducting Nb film, usually interpreted as field-induced variations of 
the critical temperature, were experimentally studied by Pannetier et al. [S^i Otani et al. [119] . Geoffroy et al. 
[120j . The oscillatory dependence of the resistivity p on the perpendicularly oriented external field Hext with a 
period AH^xt close to $o / was observed at T < Tco only when the dots had been magnetized before cooling 
(here a is the period of the magnetic dot array). The appearance of minima in the p{Hext) dependence, which 
are reminiscent of the Little-Parks oscillations for multiply-connected superconductors and superconducting 
networks (see the monograph of Berger and Rubinstein [IS] and references therein), were attributed to the 
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Array of Gd-Co dots on top of Nb film 




Figure 9. (color online) Phase transition line of a superconducting Nb film with a ferromagnetic GdssCogT 
particle array (the period 4/im), adapted from Otani et al. 1119| . 

variation of the critical temperature, Tc Tc{Hext) due to the fluxoid quantization in each magnetic unit cell. 
Tunable field-induced superconductivity for arrays of out-of-plane magnetized dots 

The fabrication of periodic arrays of multilayered Co/Pd and Co/Pt magnetic dots with out-of-plane magnetic 
moments allows one to observe both the matching phenomena and the modified OP nucleation due to the 
field-compensation effect. The stray field induced by positively magnetized dots has a positive z— component 
of the magnetic field under the dots and a negative one in the area in between the dots. Thus, the magnetic 
field in the region in between the dots will be effectively compensated at Hext > 0, stimulating the appearance 
of superconductivity at nonzero Hext values (magnetic field-induced superconductivity). Lange et al. |121j 
demonstrated that the Tc maximum is located at H^xt — for the demagnetized magnetic dots and it is shifted 
towards a certain Hext.n which depends on the dot's magnetization [see the panel (a) in Fig. [10]. This quantized 
shift of the Tc was attributed to the field compensation in the interdot areas accompanied by an annihilation 
of the interstitial antivortices under the action of the external field, since (i) the number of antivortices is 
determined by the magnetic moment of the dots; (ii) the interstitial antivortices can be fully cancelled only at 
the matching fields Hcxt.n — n^o/a^ (Milosevic and Peeters j l22j ). Thus, the appearance of periodic kinks in 
the Tc{Hcxt) phase boundary with a period coinciding with the first matching field Hi can be associated with 
the fiuxoid quantization, confirming that superconductivity indeed nucleates in multiply connected regions of 
the film. 

The results of further investigations on similar hybrid systems (an array of micron-sized Co/Pt dots on top 
of an Al film) was presented by Gillijns et al. [1231 1124|, I125j . As a consequence of the rather large diameter 
of the dots, the demagnetized dot's state microscopically corresponds to a magnetic multidomain state with 
very weak stray field. As it was demonstrated in Ref. |123| . the remanent magnetization of the dots, which 
were initially demagnetized, depends monotonously on the maximal applied field (excursion field) Hret- Thus 
the total remanent magnetic moment of the dot becomes variable and tunable, hereby changing the influence 
of the ferromagnet on the superconductor in a continuous way. It was found that a gradual increase of the 
dot's magnetization from zero to a certain saturated value results (i) in a quantized displacement of the main 
Tc maximum towards nHi {n is integer) due to the quantized character of the field-induced superconductivity 
[the panel (b) in Fig. [10]; (ii) in an enhancement of the local Tc maxima, attributed to the formation of 
a commensurate vortex phase at discrete matching fields, which becomes more pronounced as compared to 
Ref. [T2T] . 

The effect of changing the average remanent magnetization Mrem and the radius Rf oi the magnetic dots 
on the superconducting properties of an Al film deposited on top of a periodic array of such dots was studied 
by Gillijns et al. |124[ 1125] . Indeed, once the dot's magnetization becomes saturated, the only way to further 



CONTENTS 



21 



Array of Co/Pd dots on top of Pb film Array of Co/Pt dots covered by Al film 
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Figure 10. (color online) Field-induced superconductivity in a superconducting film with an array of magnetic 
dots: 

(a) The Tc{Hext) dependences obtained for a Pb film after demagnetization of the ferromagnetic dot array (the 
brown central curve marked by circles), saturation of the dots in a large positive Hext (the red right curve 
marked by diamonds), and saturation in a large negative Hext (the blue left curve marked by squares), adapted 
from Lange et al. |121| . The period of the lattice was 1.5 ^tm. The arrows depict the corresponding matching 
fields. 

(b) Superconducting transition Tc{H) of an Al film for different magnetic states of the square array of the 
ferromagnetic dots of the period 2 /xm, adapted from Gillijns et al. |123l . By increasing the magnetization a 
clear shift of Tc{Hext) and a decrease of T™"^ is observed. 



increase the magnetic fiux from each magnet can be achieved by increasing the lateral size of the dots. It was 
experimentally found that the larger the Rf value, the smaller the necessary AM needed to shift the main 
Tc{Hext) maximum by one matching field. 

Both the field compensation and matching effects in plain Al films with a square array of ferromagnetic 
Co/Pt disks were investigated by Gillijns et al. [HS] and Aladyshkin et al. |107j . Due to the presence of the 
out-of-plane magnetized dots, there are three different areas, where the OP can potentially nucleate: above the 
positive or negative domains, inside the magnetic dot, and in between the dots. In the demagnetized state 
the interdot field is close to zero, therefore superconductivity starts to nucleate at this position at relatively 
low magnetic fields, resulting in an almost linear phase boundary centered at He^t — 0. By magnetizing the 
dots positively (i) the amplitude of the field between the dots increases negatively and (ii) the typical width 
of the positive domains becomes larger than that for negative domains. Therefore the peak, associated with 
the OP localization between the dots, shifts towards positive fields. In addition, a second local Tc maximum, 
corresponding to the appearance of superconductivity above the broader positive domains, appears, while 
the OP nucleation above narrower negative domains is still suppressed. For negatively magnetized dots the 
reversed effect occurs. It is important to note, that the amplitude of the main Tc peak, corresponding to the 
OP nucleation between the dots, remains almost constant, since the mentioned area of the OP localization is 
almost independent of the dot's magnetic state (Fig. [TT|) . 

Individual ferromagnetic dots above/inside superconducting films 

Marmorkos et al. [126] studied a possibility to create giant vortices by a ferromagnetic disk with out-of-plane 
magnetization embedded in thin superconducting film within full nonlinear, self-consistent Ginzburg-Landau 
equations. Later using the same model Milosevic and Peeters [1271 11281 I129j considered the formation of 
vortex-antivortex structures in plain superconducting films, infinite in the lateral direction, in the field of an 
isolated ferromagnetic disk with out-of-plane magnetization within the full nonlinear Ginzburg-Landau theory. 
Antivortices were shown to be stabilized in shells around a central core of vortices (or a giant vortex) with 
magnetization-controlled "magic numbers" (Fig. I12p . The transition between the different vortex phases while 
varying the parameters of the ferromagnetic dot (namely, the radius and the magnetization) occurs through the 
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Figure 11. (color online) The phase boundaries Tc^H^^^t) for an S/F hybrid, consisting of an Al film and 
an array of magnetic dots, in the demagnetized state, in the completely magnetized state in positive direction 
as well as in several intermediate magnetic states, adapted from Gillijns et al. |96l and Aladyshkin et al. 
11071 . The period AHext of the Tc oscillations, which are distinctly seen in the curve corresponding to the 
magnetized states, is equal to 5.1 Oe and it exactly coincides with the matching field, i.e., AHext = ^o/'S, 
where $o = 2 ■ 10~^ Oe-cm-^ is the flux quantum and 5 = 4 fim^ is the area of the unit cell. Note that field and 
temperature intervals shown here are much broader than those presented in Fig. [To] (b). 



creation of a vortex-antivortex pair under the magnetic disk edge. 
2.6. Mesoscopic S/F hybrids: theory and experiments 

In all the description of nuclcation of superconductivity so far, we have ignored the effects of the sample's 
borders. It is well known that the OP patterns in mesoscopic superconducting samples with lateral dimensions 
comparable to the superconducting coherence length and magnetic penetration depth, is substantially influenced 
by the geometry of the superconductor (see review of Chibotaru et al. [40] and references therein). As a 
result, the presence of the sample's boundaries allows the appearance of exotic states (giant vortex states, 
vortex clusters, shell configurations etc.), otherwise forbidden for bulk superconductors and non-patterned plain 




Figure 12. (color online) Contour-plots of the Cooper pair density IV'Pi illustrating the appearance of vortex- 
antivortex shell structures in large-area superconducting film in the field of perpendicularly magnetized disk 
for different magnetic moments m of the ferromagnetic particle: raj rat) = 25 (a), 29 (b), 35 (c) and 38 (d), 
where rrto = -ffc2 5'^i by courtesy of M.V. Milosevic. The highest values are shown in lighter shades and 
the lowest densities in darker shades. It is important to note that only the central part of the superconducting 
film (45 5 X 45^) is shown here, black solid line schematically depicts the edge of the magnetic dot. Red circles 
correspond to the vortex cores, while blue squares mark the position of antivortices. The number of vortices iV„ 
and antivortices Nav are indicated on the plots. The simulations were performed for the following parameters: 
K = A/5 = 1.2, the lateral size of superconducting sample is 256 256 5, the radius of ferromagnetic disk is 
4.53 5, and the thicknesses of superconductor and ferromagnet are equal to 0.1 ^. 
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superconducting films. Since, as we have pointed out above, a nonuniform magnetic field is an alternative way 
to confine the superconducting OP in a certain H^xt and T range, mesoscopic S/F hybrids seem to be of interest 
for studying the interplay between different mechanisms of confinement of the superconducting condensate. 

It is important to note, that the screening effects can still be omitted provided that the lateral size of the 
thin superconducting sample is smaller than the effective penetration depth = ^^/Ds- In this case the 
self-interaction of the superconducting condensate can be taken into account solving the nonlinear decoupled 
GL equation 

n 2 



^ 27rz . , , 



^_^+|^|2^ = 0, (16) 



where the vector potential distribution is given by the external sources and the ferromagnet only [see Eq. ^] . 
Interplay between different regimes of the OP nucleation 

As we anticipated above, a very interesting phenomenon in mesoscopic S/F hybrids is the interplay between 
competing regimes of the OP nucleation, which can be clearly seen in the case of a small-sized magnetic 
dot of radius Rj placed above a mesoscopic superconducting disk, i?/ <C i?s- Indeed, the \tp\ maximum can be 
generally located either at the central part of the superconducting disk, close to the magnetic dot (magnetic-dot- 
assisted superconductivity) or at the outer perimeter of the superconducting disk (surface superconductivity). 
For a positively magnetized dot the regime of the magnetic-dot-assisted superconductivity, associated with 
the appearance of superconductivity in the region with compensated magnetic field, can be realized only for 
Hext < 0. At the compensation field (ji/ea:*! — -Bo) and provided that ^yW(J{2^^Bo) <^ Rg the enhancement of 
the z— component of the field near the disk edge acts as a magnetic barrier for the superconducting condensate 
and it prevents the edge nucleation of superconductivity even in small-sized superconductors {Bq being the 
maximum of the self-field of the magnetized dot). The OP nucleation near the magnetic dot becomes possible, 
if the critical temperature Tc^"^ of the formation of localized superconductivity with the OP maximum at the 
superconducting disk center 

2^0) 27r^2 

exceeds the the critical temperature for the edge nucleation regime 
1- ^ ^0.59^1 1^^,, I, 

corresponding to the critical field of surface superconductivity Hc3 = 1.69 -ffc2 [1311311133]. Due to the different 
slopes dTc°^ /dHext and dTcs/dHext and the offsets, one can conclude that the edge OP nucleation regime 
apparently dominates both for positive and large negative Hgxt values. Only in the intermediate field range the 
highest critical temperature corresponds to the formation of superconductivity near the magnetic particle. 

Little-Parks oscillations in mesoscopic samples 

The nucleation of superconductivity in axially-symmetrical mesoscopic S/F structures (e.g., superconducting 
disks or rings in the field of a perpendicularly magnetized ferromagnetic circular dot) were studied theoretically 
by Aladyshkin et al. [53] , Cheng and Fertig [l30j , Milosevic et al. |1311 1132] , and experimentally by Golubovic 
et al. [133[ I134[ 11351 1136j and Schildermans et al. [137] . Due to the cylindrical symmetry of the problem, 
superconductivity was found to appear only in the form of giant vortices 'ip{r, 9) — fhir) exjp{iL0), where L is 
the angular momentum L of the Cooper-pairs (vorticity) . The appearance of vortex-antivortex configurations 
in superconducting disks of finite radius at temperatures close to is possible, although these states were 
predicted to be metastable states. 

The observed periodic cusp- like behavior of the Tc{Hext) dependence was attributed to the field- induced 
transition between states with different vorticity similar to that of mesoscopic superconductors in a uniform 
magnetic field. However, the stray field, induced by the magnetized dot, was shown to be responsible for a 
peculiar asymmetry of the oscillatory Tc{Hext) phase boundary and a shift of the main Tc maximum towards 
nonzero Hext values [531 11331 11341 1135j . The mentioned abrupt modification of the preferable nucleation regime 
when sweeping Hf-xt can lead to a double change in the slope of the Tc{Hext) envelope from Tco/h'^ to 
Tco/H^f [Ml [137] (see Fig. The restoration of the slope close to T^o/Hf^ can be interpreted as an effective 
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Figure 13. (color online) The phase transition lines Tc{H^xt), obtained experimentally for the mesoscopic 
S/F hybrid system with a 0.1 -Rn criterion for three different magnetic states (completely magnetized, partly 
magnetized and demagnetized states), i?„ and Hm being the normal-state resistance and the magnetizing field, 
adapted from Schildermans et al. 11371 . The considered S /F system consists of superconducting Al disk covered 
by ferromagnetic Co/Pt multilayered film of the same radius Rs = Rf = 0.825 ^m. Black solid line represents 
the Tc{Hext) dependence for the reference Al disk of the same lateral size. 



elimination of the boundary effects in mesoscopic S/F samples at the compensation field (near the main Tc 
maximum). Interestingly, the nonuniform magnetic field can be used to control the shift in the field dependence 
of the maximal critical current Ic{Hext) for a bias current flowing through the superconducting loop, which 
allows one to tune the internal phase shift in superconducting networks (Golubovic et al. [136] V 

It is important to note that the periodicity AH^xt of the Little-Parks oscillations in the Tc{Hext) dependence 
is explicitly given by the area where the superconducting OP is confined. For edge nucleation only the area 
enclosed by the superconductor determines the period of the oscillations, which can be roughly estimated 
as AHext — ^o/R^- However in the case of the magnetic-dot- assisted nucleation (in the vicinity of the 
compensation field) the area of the OP localization is determined by the spatial characteristics of the nonuniform 
magnetic field (either the dot' radius Rf or the vertical separation between the dot and the superconductor 
Zd), therefore AH^xt — 'I'o/ niaxji?^, Zj}. As a consequence, the change of the nucleation regimes manifests 
itself as an abrupt modification of the oscillatory Tc{Hext) dependence. In particular, both the amplitude and 
the period of the Little-Parks oscillations become much larger, provided that (Zd^Rj) <C Rs and the OP wave 
function localizes far from the sample's edges (Aladyshkin et al. [53], Carballeira et al. [138| ). 

Symmetry-induced vortex-antivortex patterns 

Mesoscopic S/F hybrids of a reduced symmetry (e.g., structures consisting of a superconductor/ferromagnet 
disks and/or regular polygons) represent nice model systems for studying symmetry- induced phenomena. The 
formation of different vortex-antivortex configurations was studied theoretically by Carballeira et al. |138j 
and Chen et al. [139) for mesoscopic superconducting squares with a circular ferromagnetic dot magnetized 
perpendicularly, and experimentally by Golubovic et al. |140) for a superconducting Al disk with a magnetic 
triangle of Co/Pt on top. It was shown that the symmetry-consistent solutions of the Ginzburg-Landau 
equation^, earlier predicted for mesoscopic superconducting polygons by Chibotaru et al. [ID], are preserved 
for regular superconducting polygons in the stray field of ferromagnetic disk. However, since spontaneously 
formed vortices and antivortices interact with the magnetic dot in a different way, it leads to a modification 
of the symmetry-induced vortex patterns (see Fig. I14p . In particular, the dot can be used to enlarge these 
vortex-antivortex patterns, thus facilitating their experimental observation with local vortex- imaging techniques 

t By symmetry-consistent solutions we mean those vortex configurations reflecting the symmetry of the problem. For instance, in 
a mesoscopic superconducting square with vorticity L = 3, the state consisting of four vortices and a central antivortex may have 
lower energy than the configuration of three equidistant vortices, which breaks the square symmetry (Chibotaru et al. |141| ). 
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Figure 14. (color online) A comparison between the vortex-antivortex patterns that can be observed for L = 3 
(a) with and (b) without the magnetic dot on top of the superconducting square, a is the lateral size of the 
sample, adapted from Carballeira et al. 11381 . As can be seen, the vortex-antivortex pattern rotates 45° and 
expands dramatically in the presence of the magnetic dot. 

( "magnetic lensing" ) . 
Embedded magnetic particles 

Doria |1421 1143j and Doria et al. [1441 1145] studied theoretically the formation of vortex patterns induced 
by magnetic inclusions embedded in a superconducting material but electrically insulated from the multiply- 
connected superconductor. Since in the absence of an external field, flux lines should be closed, vortex lines 
are expected to start and end at the magnetic inclusions. The calculations, performed in the framework of GL 
theory for a mesoscopic superconducting sphere with a single magnetic point like particle in its center, reveal 
that the confined vortex loops arise in triplets from the normal core when the magnetic moment reaches the 
scale defined by mo = ^o^/2Tr. Therefore a vortex pattern is made of confined vortices, loops and also broken 
loops that spring to the surface in the form of pairs. This vortex state provides a spontaneous vortex phase 
scenario for bulk superconductors with magnetic inclusions, where the growth of the vortex loops interconnects 
neighboring magnetic inclusions. 

Effect of the finite thickness of the superconducting film 

The effect of the finite thickness of the superconducting film on the Tc{Hext) dependence was studied 
theoretically by Aladyshkin et al. [M], Aladyshkin and Moshchalkov |95j, and Schildermans et al. |137j . 
It follows from Maxwell's equations that the faster the spatial variation of the magnetic field in the lateral 
(a;, J/)— directions, the faster the decay of the field in the transverse z— direction. However the OP variation 
along the sample's thickness can be effectively suppressed by a requirement of vanishing the normal derivative 
dip / dn of the superconducting OP at the top and bottom surfaces of the superconductor. Indeed, the 




Figure 15. (color online) Confined vortex loops arise in sets of threes inside a mesoscopic superconducting 
sphere (radius 15 J), as shown here for two consecutive values of the point-like magnetic moment that occupies 
its center: (a) m/mo = 15; (b) m/mo = 20, by the courtesy of M. M. Doria. 
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external field applied along the z— direction makes the OP variations over the superconducting film thickness 
more energetically unfavorable than in the lateral direction, especially for H^xt values of the order of the 
compensation field or higher. The field dependence Tc{Hext) for rather thick superconducting films in the 
high-field limit is similar to the phase transition line described by Eq. ([7]). This behavior results from the 
effective averaging of the inhomogeneous magnetic field by the quasi-uniform OP wave function over the sample 
thickness, which substantially weakens the effect arising from the field modulation in the lateral direction. 
In other words, only the lateral inhomogeneity leads to the anomalous Tc{Hf,xt) dependence and reentrant 
superconductivity in particular, while the transverse field non-uniformity masks this effect. Interestingly, the 
location where superconductivity starts to nucleate shifts towards the point with zero total magnetic field (not 
only z— component of the magnetic field) as the thickness of superconductor increases (Aladyshkin et al. [54]). 
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3. Vortex matter in non-uniform magnetic fields at low temperatures 

3.1. London description of a magnetically coupled S/F hybrid system 

In section [2] we have reviewed the theoretical modelling within the GL formalism and the experimental data 
concerning the superconducting properties of S/F hybrids rather close to the phase transition line Tc{Hf.xt)- In 
this section we consider the properties of S/F hybrids for a fully developed superconducting OP wave function, 
i.e. for T <^ Tcq. In this limit, the superconducting properties of the S/F hybrids can be correctly described by 
the London theory, omitting any spatial variations of the OP. As before, we assume that the coercive field of the 
ferromagnet is much higher than the upper critical field of the superconductor. This guarantees that, during 
the investigation of the superconducting properties, no changes in the ferromagnetic element(s) will occur. 

At rather low temperatures one should take into account that the magnetic field, induced by screening 
(Meissner) currents or by vortices, can no longer be neglected and will strongly interact with the ferromagnet. 
The free energy functional of the S /F hybrid can be written in the following form [44[ I45j 

Gsf = Gso + G„, + J |^(rotB)' + ^-B.M-^:^| dV. (17) 

V 

As before, Gso is the self-energy of superconductor, while the term Gm depends on the magnetization of 
ferromagnet only. Assuming that Gm is constant for a given distribution of magnetization and minimizing this 
functional Ggf with respect to the vector potential A, one can obtain the London-Maxwell equation: 

B + X^ rot rot B = $o X! ~ + '^^^'^ ™^ ™* ^' 

i 

where summation should be done over all vortices, positioned at points {Ri,.i} and 5{r) is the Dirac delta- 
functiorjjl. In Eq. psp we assumed that each vortex line, parallel to the z— axis and carrying one flux quantum, 
generates a phase distribution 8 = with rot(VO) = rot ([zq x ro]/r) — 6{r) Zq [here {r,ip,z) are cylindrical 
coordinates with the origin chosen at the vortex] . It should be noticed that the solution of Eq. (fTS]) gives the 
magnetic field distribution for a given vortex configuration, and in order to find the field pattern one should 
find the minimum of the total free energy Gsf with respect to the vortex positions {Rt,^i}. 

3.2. Interaction of a point magnetic dipole with a superconductor 

Next we consider the generic problem of the interaction of a superconductor with a point magnetic dipole 
positioned at a height Zd above the superconducting sample. Introducing the OP phase, which determines the 
density of supercurrents js = c/(47rA^) [$o/(27r)V8 — A], the free energy functional Eq. (|17p at Hext = can 
be rewritten as Gsf = Gso + G„i + Gg , where 

- / { sJa^^^^ • - le^A^^^^ • A) - i B . m} (19) 

similar to that obtained by Erdin et al. |146| . Due to the linearity of Eq. we can introduce the field 

B^ — rot A„ , generated by a single vortex line positioned at r = , and the field B„j — rot A,„ , corresponding 
to the screening (Meissner) current distribution induced by the magnetic dipoleQ. As a result, for the particular 
case of a point magnetic dipole, M = mo (5(r— R^), the energy of the system Gs depending on the supercurrents 
can be represented as a sum of the self-energy of the vortex el"-' — ($o/47rA)^ InA/^, which is independent on 
the dipolar moment, and a term Gint, responsible for the interaction between the dipole and superconductor 
(e.g., Wei et al. [147], Carneiro 148J): 

G„t = - i mo • B„ (Rd) - mo • B„ (R^) . (20) 

— oo — oo 

t We follow the standard definition of tiic 5— function: 5{x, y) = for any x =^ 0, y =^ 0, and J J 5{x, y) dxdy = 1. 

— oo — oo 

X Tliese fields Bu and Bm are tiie solutions of tiie following differential equations 
Bt, + rot rot Bj, = <I>o <5(r - R„) zq, 
Bm + rot rot Bm = iirX^ rot rot M. 
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The first term in Eq. ([^0]) corresponds to the interaction between the magnetic dipole and the local field at the 
dipole's position Bm(R(i) due to the screening currents. The second term in Eq. (PO)) gives the energy of the 
interaction between the vortex and the magnetic dipole, which consists of two parts: (i) the "hydrodynamical" 
interaction between circulating supercurrents due to the vortex on one hand and the dipole on the other hand; (ii) 
the interaction between the stray field of the vortex with the magnetic moment. Interestingly both contributions 
turn out to be equal to (—1/2) mg • Bu(Rc;). 

Interaction between a point magnetic dipole and the Meissner currents 

In the case of a superconductor without vortices [the so-called Meissner state, rot(V0) = 0], the interaction 
between the dipole and the superconductor is given by Um = (^1/2) kiq • B,„(Rrf). It can easily be seen that 
this energy contribution is always positive, since the screening currents induced by a dipole generate a magnetic 
field which is opposite to the orientation of the dipole. This is in principle just a reformulation of the fact that 
a magnetic dipole will be repelled by a superconducting film with a force f = —dUm/dRd- The related problem 
of levitating ferromagnetic particle over various superconducting structures was considered theoretically by Wei 
et al. [HT], Xu et al. [Hn], Coffey [Unillll] and Haley [Til[T13]. In particular, in the limiting case Z^/A > 1 
the magnetic levitation force was found to vary linearly with A, regardless of the shape of the magnet. The 
fact that A (as well as the vertical component of the force) shows an exponential temperature dependence for 
s— wave superconductors, linear— T for d— wave superconductors, and quadratic— dependence in a wide low- 
temperature range for materials with s + id symmetry of the gap, can assist in distinguishing the type of pairing 
in real samples based on magnetic force microscopy (MEM) measurements [149] . 

Creation of vortices by the stray field of a magnetic dipole 

In reality, magnetic dipoles which have a high dipolar moment or are localized rather close to the superconductor 
can suppress the Meissner state due to the dipole's stray field. The possible appearance of a vortex would change 
the total energy of the S/F hybrid system by Gs = + Smv, where emv oc —mo accounts for the dipole- vortex 
interaction. Clearly, one can make Gg < by increasing toq and therefore it will become energetically favorable 
to generate vortices in the system. Such a scenario for the formation of the mixed state was considered, e.g., 
by Wei et al. [1471 1154] for a vertically magnetized dipole placed above a thin superconducting film. As a 
consequence of Eq. (|20p . the appearance of a vortex line in the superconductor drastically changes the resulting 
force acting on the magnetic dipole (Xu et al. [149j . Wei et al. |147[ 1154] ). 

It is important to notice that in the case of superconducting films infinite in the lateral direction and 
cooled in zero magnetic field, a spontaneously induced vortex should be accompanied by an antivortex in order 
to provide flux conservation imposed by Maxwell's equation^. The destruction of the Meissner state in a 
zero-field-cooled superconducting film was theoretically studied by Melnikov et al. [155] and Aladyshkin et 
al. [156j . In this case the formation of the mixed state occurs via the penetration of vortex semi- loops which 
split into vortex-antivortex pairs. The local suppression of the Bean-Livingston energy barrier |157| . which 
controls the process of vortex penetration, takes place when the screening current density will be of the order 
of the depairing current density. The threshold distance = Z^iT), corresponding to the suppression of 
the energy barrier can be experimentally detected by measuring a nonzero remanent magnetization as long 
as pinning is relevant. Interestingly, such non-contact technique allows one to estimate the depairing current 
density and its temperature dependence in thin superconducting YBa2Cu307 films jl55| 11561 1158j . Since the 
surface energy barrier in superconductors with a flat surface is known to be suppressed by applying an external 
field on the order of the critical thcrmodynamical flcld Hem = *i'o/(2\/27rA^) [43^, one can get a rough criterion 
for the persistence of the vortex- free state in the presence of a magnetic dipole: mo/Z^ < Hem- A different 
criterion should be obtained in case the dipole is already present close to the superconductor and the whole 
system is cooled down below the superconducting temperature. Under this field-cooled condition a lower critical 
magnetization to induce a vortex-antivortex pair is expected. 

The problem of the formation of vortex-antivortex pairs in superconducting films in the presence of vertically 
and horizontally magnetized dipole at H^xt = was considered theoretically by Milosevic et al. |159[ I160j , 
Carneiro ,148j. It was shown that an equilibrium vortex pattern could consist of spatially separated vortices 
and antivortices (for rather thin superconducting films) or, when the superconducting film thickness increases 

t In mesoscopic superconducting systems, the returning fiux lines generated by an out-of-plane magnetized dipole can bypass the 
border of the sample and the existence of an antivortex is not required. 
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and becomes comparable to the London penetration depth, curved vortex Unes, which start and terminate at 
the surface of the superconductor. 

Magnetic pinning 

Now we will discuss the properties of S /F hybrids where the magnetization of the ferromagnet and the distance 
between ferromagnet and superconductor are assumed to be fixed resulting in a constancy of the interaction 
energy with the screening currents induced by the ferromagnet. In this case, any variation of the free energy 
of the S/F hybrid in the presence of vortices (either induced by the ferromagnetic element or by external 
sources) can be attributed to the rearrangements of the vortex pattern. The part of the interaction energy Gint, 
proportional to the magnetization of the ferromagnet and sensitive to the vortex' positions, is usually called the 
magnetic pinning energy Up. 

In order to illustrate the angular dependence of the interaction between a point magnetic dipole of fixed 
magnetization and a vortex, we refer to the following expression (Carneiro jl48j, Milosevic et al. \159\ 1160] ) 



$2 

^.K,0) = ^ 



' ' ' I — — — cos ( 



(21) 



obtained for a thin superconducting film Ds <^ X and for vortex-dipole separations smaller than the effective 
penetration depth X2D = X^/Dg. Here mo^^ and too,z are the in-plane and out-of-plane components of the 
dipolar moment, r± ~ \/ + y^, R-d = (0, 0, Zd) is the position of the dipole, is the angle between the a:-axis 
and the vector position of the vortex in the plane of the superconductor r^. The resulting pinning potentials 
for an in-plane and out-of-plane magnetized dipole, derived from Eq. (|2ip . are shown in Fig. 1161 The pinning of 
vortices (and antivortices) in the superconducting films of a finite thickness on the magnetic dipole at iJext = 
was analyzed by Milosevic et al. |1591 I160j and Carneiro [1481 1161] . In addition, Carneiro considered the 
magnetic pinning for the case of a dipole able to rotate freely in the presence of the external magnetic field 
[m UMl HM] or external current [T^ . 

In Ref. |148j it was shown that this magnetic pinning potential has a depth of toq / (47rA) , penetrates a 
distance A into the film whereas its range parallel to the film surface is a few times A. This finding points out 
the relevance of the penetration depth A to characterize the purely magnetic pinning potential. Since typical 
A values are similar for a broad spectrum of superconducting materials, this suggests that magnetic pinning 
represents a promising way of increasing the critical current not only in conventional superconductors but also 
in high- Tc superconductors. Notice that the fact that the magnetic pinning range is determined by A sets a limit 
for the minimum distance between magnetic particles, beyond which the vortex lines cannot resolve the field 
modulation and therefore the pinning efficiency decreases. Having in mind that many practical applications 




Figure 16. (color online) The spatial dependence of the energy of a single vortex in a field of point magnetic 
dipole, M = mo<5(a::, y, z — Z^), calculated according to Eq. 1211 1 for the height Z^jX = 1 and the dipole strength 
mo/(*oA) = 10: 

(a) vertically magnetized dipole, mo = niQZo, 

(b) horizontally magnetized dipole, mo = moXQ. 
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typically involve high-Tc superconductors with ^ ^ A, it becomes clear that this maximum density of pinning 
sites to trap individual vortices is much lower than that limited due to core pinning. 

It is important to mention that real nanoengineered ferromagnetic elements are quite far from point dipoles 
but rather consist of extended volumes of a magnetic material, such as magnetic dots and stripes. In this case, 
due to the principle of superposition, the magnetic pinning energy needs to be integrated over the volume of 
the ferromagnet Vf to determine the interaction between a vortex line and a finite size permanent magnet 



The magnetic pinning in an infinite superconducting film produced by individual ferromagnetic objects 
of finite size was considered for the following cases: a ferromagnetic sphere magnetized out-of-plane (Tokman 
[165j ). magnetic disks, rings, rectangles, and triangles magnetized out-of-plane (Milosevic et al. |160[ I166| l. 
magnetic bars and rectangular dots magnetized in-plane (Milosevic et al. [167] ). circular dots magnetized either 
in-plane and out-of-plane (Erdin et al. |146| . Erdin |168| ). circular and elliptic dots and rings magnetized 
out-of-plane (Kayali [1691 1170j . Helseth [171j ), columnar ferromagnetic rod (Kayali |172| 1. The analysis of 
the pinning properties of a vortex in a semi-infinite superconducting film due to magnetic dots was done by 
Erdin (173j . Similar to the case of a point magnetic dipole, the extended magnetic objects are able to generate 
vortex-antivortex pairs provided that the size and the magnetization are large enough. 

We would like to note that the attraction of a superconducting vortex to a source of inhomogcneous 
magnetic field (e.g., coil on a superconducting quantum interference device or magnetized tip of a magnetic force 
microscope) makes it possible to precisely manipulate the position of an individual vortex. Such experiments 
were performed by Moser et al. |174| . Gardner et al. |175| . Auslaender et al. |176j for high— Tc superconducting 
thin films and single crystals in intermediate temperatures when intrinsic pinning become weaker. This technique 
apparently opens unprecedent opportunities, e.g., for a direct measuring of the interaction of a moving vortex 
with the local disorder potential and for preparing exotic vortex states like entangled vortices (Reichhardt [177] ). 

Equation (|22p shows that the pinning potential does not only depend on the size and the shape of the 
ferromagnetic elements but also on the particular distribution of the magnetization (i.e. on their exact magnetic 
state). It is expected that the average pinning energy is less efficient for magnetic dots in a multidomain state 
whereas it should reach a maximum for single domain structures. In other words, if the size of the magnetic 
domains is small in comparison with A or the separation distance Zd^ then a vortex line would feel the average 
field emanating from the domains and the magnetic pinning should be strongly suppressed. This flexibility of 
magnetic pinning centers makes it possible to tune the effective pinning strength, as we shall discuss below. 

The question now arises whether the pinning potential produced by a magnetic dipole will remain efficient 
when a vortex-antivortex pair is induced by the magnetic dipole. In order to answer this question it is necessary 
to minimize the mutual interaction between the induced vortex-antivortex pair, the magnetic element and the 
test vortex generated by the external field. The magnetic moment - test vortex interaction is a linear function 
of Too always favoring the test vortex to sit on the positive pole of the magnet, whereas the induced currents - 
test vortex interaction will follow the Little-Parks oscillations (due to the creation of vortices by the dipole) . For 
the case of an in-plane dipole, Milosevic et al. [167j demonstrated that all these terms produce a subtle balance 
of forces which lead to a switching of the optimum pinning site from the positive to the negative magnetic pole, 
as Too is increased. For the case of out-of-plane magnetic dipole, once a vortex-antivortex pair is induced, the 
test vortex will annihilate the antivortex leaving a single vortex on top of the magnetic dot (Gillijns et al. [123] 1. 

3.3. Magnetic dots in the vicinity of a plain superconducting film 

The rich variety of possibilities considered theoretically in the previous section for an individual ferromagnetic 
element in close proximity to a superconducting film represents an experimental challenge in part due to the 
difficulties associated with recording small induced signals. A very successful way to overcome this limitation 
consists of studying the average effect of a periodic array of dots at the expense of blurring sharp transitions, 
such as the vortex-antivortex generation, or inducing new collective phenomena associated with the periodicity 
of the magnetic templates. An important experimental condition that should be satisfied in order to justify 
the analogies drawn between individual dipoles and arrays of dipoles, is the lack of magnetostatic interactions 




(22) 
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Array of Co dots covered by Pb film 




Figure 17. (color online) Upper half of the magnetization loop M vs. H^^t at T/T^q = 0.97, for a 
superconducting Pb film (50 nm thickness) on top of a triangular lattice of Au/Co/Au dots (period 1.5 /im 
corresponding to a first matching field of 9.6 Oc) before and after magnetizing the dots (filled and open symbols, 
respectively), and for a reference Pb 50 film (solid line), adapted from van Bael et al. |196| . The curve for the 
magnetized dots is slightly shifted upwards for clarity. 

between neighboring dipoles. In other words, it is necessary to ensure that the field generated by a dipole at the 
position of its nearest neighbors hes below the coercive field of the chosen ferromagnetic materials (Cowburn et 
al. [T7i[179], Novosad et al. [1801 fTST] ). 

Commensurability effects in S/F hybrids with periodic arrays of magnetic dots 

In the early 1970's, Antler [1821 1183] proposed that a periodic array of ferromagnetic particles should give rise 
to an enhancement of the critical current of the superconducting material. Recent developments on lithographic 
techniques have made it possible to prepare superconducting structures containing a regular array of magnetic 
dots (Co, Ni, Fe) at submicrometer scale of desirable symmetry in a controlled way (Otani et al. [119] . Geoffroy 
et al. [110], Martin et al. [HI (113 HHS dHT] , Morgan and Ketterson [UHldHn], Hoffmann et al. [THD], Jaccard 
et al. [HI], Villegas et al. [Ml UMl HM], Stoll et al. [US], van Bael et al. [Ml (Ml [Ml [Ml 1^ - 
van Look et al. |202| ). The S/F hybrids containing periodic arrays of magnetic elements with out-of-plane 
magnetization (multilayered Co/Ft and Co/Pd structures) were fabricated and investigated by van Bael et al, 
[Ml [Ml UnOl Uni 1203 1201] and Lange et al. [103 HOI]. In all these works, it was found that the presence 
of the lattice of magnetic dots leads either to (i) a resonant change in the magnetoresistance p{Hf.^t) and the 
appearance of pronounced equidistant minima of resistivity with the period Hi — $o / Sq determined by the 
size of the magnetic unit cell So [m [1201 HM [113 UMl [Ml [113 [Ml [Ml [Ml [Ml [Ml [M]; (u) or to the 
presence of distinct features as peaks or plateaus in the field dependence of the critical current Ic{Hext) or in 
the magnetization curve M{Hext ) [Ml [Ml [Ml [Ml [Ml [Ml [103 HOIl HOI HOI [103 [103 [100] ■ 

Such matching effects are also very well known for spatially-modulated superconducting systems with 
antidot lattices without ferromagnetic constituents (e.g., Baert et al. [51], Moshchalkov et al. [53], Hebard et al. 
|207| . Rosseel et al. [2081 , Harada et al. [209], Metlushko et al. |210j ). These periodic anomalies are commonly 
explained in terms of commensurability effects between the vortex lattice governed by the external field and the 
artificially-introduced pinning potential at an integer number of vortices per unit cell at H^^t — ±nHi (with n 
integer). Typically these lithographically defined arrays are limited to a minimal period of the unit cell of the 
order of a few hundred nanometers, giving rise to a maximum matching field Hi ~ 10 — 10^ Oe. Alternative 
methods for introducing more closely packed particles and thus larger Hi values can be achieved by Bitter 
decoration (Fasano et al. |211[ I212j , Fasano and Menghini |213j ) or a diversity of self-assembled techniques 
(Goyal et al. [H^, Villegas et al. [215], Welp et al. [2161 [217], Vinckx et al. [218], Vanacken et al. [2T9] ). 

The field dependence of the critical current Ic{Hext) and magnetization M{Hext), can be symmetrical or 
asymmetrical with respect to H^xt = depending on the dot's net magnetization (compare Figs. [TT] and [TS]) . 
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Figure 18. (color online) Magnetization curves M vs. Hext at T = 7 K (Tco = 7.17 K, T/Tco = 0.976) for a 
superconducting Pb film (50 nm thickness) on top of a Co/Pt dot array (the period is 1.0 /xm, the first matching 
field is 20.68 Oe) with all dots aligned in the positive (a) and negative (b) direction, adapted from van Bael et 
al. [204) . 



The latter takes place for arrays of magnetic dots with out-of-plane magnetization |1981 11991 12001 12011 12031 
12041 12051 1206j . Indeed, vertically magnetized dots with average magnetic moment (111)2 > 0, similar to point 
magnetic dipoles, produce a stronger pinning potential for vortices (at H^xt > when (m) || Hext) than for 
antivortices {Hext < 0). In contrast to that, in-planc magnetized dots are able to pin vortices and antivortices 
at the magnetic poles equally well [see Eq. (PT|) and Fig. [TB]. This explains the experimentally observed field 
polarity-dependent (asymmetric) pinning for arrays of out-of-plane magnetized particles (Fig. [TS]and Fig. I19p . 

The interaction between vortices and a periodic array of hard magnetic dots on top or underneath a plain 
superconducting film within the London approximatiorU was theoretically analyzed by Helseth [225] , Lyuksyutov 

t This issue seems to be part of a more general problem of the interaction of vortex matter with a periodic potential regardless the 
nature of the pinning in the superconducting system (see, e.g., Reichhardt et al. 12201 12211 12221 [2231 12241 and references therein). 
In this review we discuss only the results obtained for the S/F hybrids, keeping in mind that similar effects can be observed for 




Figure 19. (color online) Field dependence of the critical depinning current jc, calculated for a periodic array 
of out-of-plane magnetized dots. The values of the dot's magnetization are indicated in the plot. Jgl is the 
density of the depairing (Ginzburg-Landau) current, T/TcO = 0.9, adapted from Milosevic and Peeters |122| . 
It is worth noting that (i) jc{Hext) is asymmetric similar to that shown in Fig. 1181 and (ii) the quantized 
displacement of the jc maximum toward nonzero Hext value is sensitive to the magnetization of the dots. 
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Figure 20. (color online) (a) Field dependence of the electrical resistance p for a superconducting Nb films 
covering periodic arrays of the magnetic Ni dots with different dot's diameter (indicated in the plot), but for 
the same lattice constant 400 nm, adapted from Hoffmann et al. 11901 . The curves are shifted by factors of 10 
from each other for clarity. 

(b) Dependences p{Hext) for samples with different diameters of non- magnetic Ag dots, adapted from Hoffmann 
et al. [Tool. 



and Pokrovsky [US], Sasik and Hwa [H^, Erdin [55S], Wei [UHlllSU], Chen et al. [531] • These calculations 
show that at H^xt = for out-of-plane magnetized dots, vortex-antivortex pairs can be created in thin- film 
superconductors with the vortices always sitting on top of the magnetic dot and the antivortices located in 
between the dots. For in-plane magnetized dots (or magnetic bars), the vortex and antivortex will be located 
at opposite sides of the magnetic dots as described above for individual magnetic dipoles. Unlike the case of an 
isolated dipole, the threshold magnetization value needed to create a vortex-antivortex pair is also a function 
of the period of the lattice (Milosevic and Peeters |114| ). Direct visualization of vortex-antivortex pairs via 
scanning Hall probe microscopy was achieved for a square array of in-plane dots by van Bael et al. |197| and 
for out-of-plane dots by van Bael et al. [204j and Neal et al. |232j . 

It is known that the preferred vortex configuration in a homogeneous defect-free superconducting film should 
be close to a triangular lattice because of the repulsive vortex- vortex interaction [331 ISl 131] . The artificially- 
introduced pinning appears to be the most effective provided that each vortex is trapped by a pinning center, 
i.e. when the symmetry of the pinned vortex lattice coincides with that imposed by the topology of the internal 
magnetic field. The transition between square and distorted triangular vortex lattice, induced by variation 
of the strength of the periodic pinning potential and the characteristic length scale of this interaction, was 
considered by Pogosov et al. |233| for superconductors with a square array of pinning centers. Experimentally 
the field-induced reconfiguration of the vortex lattice (from rectangular to square) for superconducting Nb films 
and rectangular arrays of circular magnetic Ni and Co dots was reported by Martin et al. |185j and StoU 
et al. |195j as an abrupt increase of the period of the oscillation in the piH^^t) dependence (resulting from 
the shrinkage of the period of the vortex lattice) and decrease of the amplitude of such oscillations (due to a 
weakening of the effective pinning) while increasing jTJea;*!- 

The dependence of the magnetic pinning in superconducting Nb films on the diameter of the Ni dots was 
studied by Hoffmann et al. [190j . They found that more minima appear in the magnetoresistance (or maxima 
in the critical current) as the lateral dot's size increases, indicating thus an enhanced pinning (the panel (a) 
in Fig. [50)1 . This effect can be caused by the two parameters which increase with the dot size: the total 
magnetic moment (m) (proportional to the dot's volume Vf = irR^jDf) and the area on the order of 7ri?y where 
superconductivity might be locally suppressed due to the high stray field or proximity effect. In addition, larger 
magnetic dots can stabilize giant vortices carrying more than one flux quantum. 

non-magnetic patterned superconductors as well. 
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Figure 21. (color online) Magnetization curves M(Hext) at T = 7.05 K (Tco = 7.20 K, T/T^o = 0.972) of a 
superconducting Pb film on top of a magnetic Co/Pt antidot lattice (the period is 1.0 fim, the first matching 
field is 20.68 Oe) after saturation in a positive field (M^ > 0, open circles) and after saturation in a negative 
field (M^ < 0, filled circles), adapted from Lange et al. [234ll235l . 



Periodic arrays of magnetic antidots 

The antipode of arrays of magnetic dots is a perforated ferromagnetic film (so-called magnetic antidots) , which 
also produces a periodic magnetic field. This system can be regarded as the limiting case of big magnetic dots 
with a diameter larger than the period of the periodic lattice. 

Magnetic antidots in multilayered Co/Pt films, characterized by an out-of-plane remanent magnetization, 
and their influence on the superconducting properties of Pb films were studied by Lange et al. [2051 12341 12351 
1236] . From magnetostatic considerations, such submicron holes in a ferromagnetic thin film generates a very 
similar field pattern as an array of magnetic dots of the same geometry, but with opposite sign. As a consequence, 
the enhanced magnetic pinning and the pronounced commensurability peaks in the M^H^xt) dependence are 
observed for the opposite polarity of the external field (i.e. at H^^t < for positively magnetized film and vice 
versa) , see Fig. [5TJ However the matching effects are considerably weakened in the demagnetized state of the 
Co/Pt multilayer with holes as compared with the demagnetized array of magnetic dots, thus indicating that 
an irregular domain structure effectively destroys a long-range periodicity [205| . 

Van Bael et al. [237 and Raedts et al. [23,8! explored perforated Co film with in-plane anisotropy. 
In this case the magnetic field distribution becomes non-trivial since such magnetic antidots effectively pin 
magnetic domain walls which generate a rather strong magnetic field. As a result, neither matching effects nor 
pronounced asymmetry were observed in the magnetization curves of the superconducting layer, but only an 
overall enhancement of the critical current after the sample was magnetized along the in-plane easy axis, in 
comparison with the demagnetized state. 

Anisotropy of the transport characteristics and guidance of vortices 

In any periodic array of pinning centra, transport properties such as magnetoresistance p{H) and the critical 
current Jc{H), exhibit a dependence not only on the absolute value of Hg^^tj but also on the direction of the 
applied transport current with respect to the principal translation vectors of the periodic pinning array (Villegas 
et al. [Ml [Ml, Soroka and Huth [239], Velez et al. [240], Silhanek et al. [241], Wordenweber et al. [242]). 
Interestingly, the direction of the Lorentz force fz, = c~^[j x B] and the drift velocity of the vortex lattice do not 
generally coincide. It was demonstrated that for rectangular arrays of magnetic dots the minimum of resistivity 
corresponds to a motion of the vortex lattice along the long side of the array cell. Such behavior was predicted 
by Reichhardt et al. [223j by numerical simulations indicating that a rectangular array of pinning centers 



CONTENTS 



35 




45 90 135 180 

p, degree 



Figure 22. (color online) Critical current, Jc, for a vortex pinned by a dipole array as a function of the angle 
(3 between the magnetic moments and the driving force, adapted from Carneiro |243| . Here Jd is the depairing 
current. 




Figure 23. (a) Magneto-optical image of a non-patterned superconducting Pb disk at T = 2 K and Hf^^t = 50 
Oe, demonstrating an isotropic flux penetration, after Gheorghe et al. I246| . White corresponds to a high local 
magnetic field and black to zero local field, (b) Magneto-optical image of a circular Pb sample decorated by 
fully magnetized Co/ Pt dots, obtained at T = 2 K and H^^t = 72 Oe, after Gheorghe et al. 12461. The external 
field is applied parallel to the dot's magnetic moment. 

(c) Magneto-optical image of fiux entry in a superconducting MoGe film at H^}^^ = 16.5 Oe, T = 4.5 K 
following the preparation of stripe domain structures in a permalloy film by turning on and off an in-plane field 
of = 1 kOe at an angle of 45° with respect to the sample edge at T > Tco, after Vlasko-Vlasov et al. 

| 247| . The brightness of the magneto-optical contrast corresponds to the vortex density. The large yellow arrow 
shows the preferential fiux entry direction coinciding with the direction of the stripe domains in the Py film. 
The thin solid line with arrow marks the sample edge, (d) Same as in panel (c) after application and switching 
off of = 1 kOe along the sample edge at T > Tco, after Vlasko-Vlasov et al. |247l . 

induces an easy direction of motion for the vortex lattice (and larger dissipation as well) along the short side of 
the array cell. 

Similar anisotropic transport properties was studied by Carneiro |243j for the case of a periodic array 
of in- plane magnetic dipoles. In order to illustrate the angular dependence of the critical depinning current 
on the angle (3 between the direction of the injected current and the magnetic moment of in-plane oriented 
dipoles we refer to Fig. [551 Interestingly, Verellen et al. [244] showed that this resulting guided vortex motion 
in square arrays of magnetic rings can be rerouted by 90° simply by changing the dipole orientation or can 
even be suppressed by inducing a flux-closure magnetic vortex state with very low stray fields in the rings. 
Similar anisotropic vortex motion was recently observed in Nb films with a periodic array of one-dimensional 
Ni lines underneath by Jaque et al. '245j. The mentioned channelling of vortices lead to an anisotropic vortex 
penetration that has been directly visualized by means of magnetooptics experiments [Gheorghe et al. [246j . 
see Fig. mb)]. 

Mechanisms of pinning in S/F hybrids 

It should be noticed that the magnetic pinning originating from the spatial modulation of the "internal" magnetic 



CONTENTS 



36 



field generally competes with so-called core pinning resulting from structural inhomogeneities in real samples 
(either regular or random defects). In addition to random intrinsic pinning, the fabrication of an array of 
magnetic particles naturally leads to an alteration of the local properties of the superconducting film (e.g., due 
to proximity effects, corrugation of the superconducting layer or local suppression of the critical temperature). 
As a consequence, both magnetic and structural modulation share the same periodicity, and a clear identification 
of the actual pinning type becomes difficult. 

A direct comparison of the pinning properties of arrays of magnetic vs. nonmagnetic dots have been 
addressed by Hoffmann et al. |190j and Jaccard et al. |191| . These reports show that even though both 
systems display commensurability features, the pinning produced by magnetic arrays of Ni dots is substantially 
stronger than that produced by non-magnetic Ag particles (Fig. I20p . In our opinion, the main issue whether the 
enhanced pinning for the sample with ferromagnetic Ni dots actually arises from purely magnetic interactions 
and not from an additional suppression of the local critical temperature, e.g. due to the enhanced magnetic 
field near magnetic dots, remains unclear. In principle, the most straightforward way to distinguish the two 
competing pinning mechanisms is the mentioned field-polarity of the magnetic pinning for the S /F hybrids with 
dots magnetized perpendicularly, i.e. exploring the broken field-polarity symmetry. 

Clear evidence of the field-polarity dependent pinning properties has been reported by Gheorghe et al. 
[246] in Pb films on top of a square array of [Co/Pt]io dots with a well defined out-of-plane magnetic moment. 
In this work the authors show that the critical current of the hybrid system can be increased by a factor of 
2 when the magnetic dots are switched from low stray field in the demagnetized state (disordered magnetic 
moment) to high stray field in the magnetized state (nearly single domain state) at temperatures as low as 
T ~ 0.3 Tco (see Fig. Additional evidence of an increase of the critical current at low temperatures (far 

from the superconducting/normal phase boundary) produced by magnetic dots was reported by Terentiev et 
al. [2481 [2491 [250] , 

Tunable pinning centers 

An apparent advantage of using magnetic pinning centra is their flexibility (tunability) in contrast to core 
pinning on structural inhomogeneities. Indeed, according to Eq. (|22p the magnetic pinning should be sensitive 
to the particular distribution of magnetization inside the ferromagnetic elements. Depending on the geometrical 
details of the dot and the magnetic anisotropy of the chosen material a huge variety of magnetic states can 
be found. For instance, domain formation is expected to be suppressed for structures with lateral dimensions 
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Figure 24. (color online) Temperature dependence of the critical current density jc, estimated from magneto- 
optical images, for a superconducting Pb film with square array of the ferromagnetic Co/Pt dots on top, 
in various magnetic states of the dots: demagnetized (o), fully magnetized parallel configuration (o), fully 
magnetized antiparallel configuration (square), partially magnetized parallel, ('m)^ = 0.2b MsVj (A) and 
(m)z = MsVf (V), adapted from Gheorghe et al. 12461 . The dashed lines are guides to the eye. 
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Figure 25. (color online) (a) Normalized magnetization M'I/Ms vs. in-plane applied field H^^^ {Ms is the 
saturated magnetization) for an array of Fe dots with average diameter of 140 nm and average interdot distance 
of 180 nm measured at T = 6 K (above the critical temperature of superconducting Al film), adapted from 
Villegas et al. |215l . Brown diamonds correspond to a virgin state of the Fe dots, the magnetic state depicted 
by red (blue) circles obtained after saturation in positive (descending branch) and negative (ascending branch) 
magnetic fields, respectively. The filled circles schematically show the regions where the magnetic vortex is 
expected to take place. Vertical dashed lines mark the coercive fields. 

(b) Normalized resistivity p vs. in-plane applied field H^^^ for the same sample at T = 1.25 K (below the critical 
temperature of superconducting Al film), p„ is the normal-state resistance. Open (blue) and filled (red) circles 
mark the curves measured from negative and positive saturation respectively, while brown diamonds correspond 
to the virgin state, adapted from Villegas et al. |215l . 



smaller than tens of nm (Raabe et al. |251j ). whereas for larger sizes the magnetic sample breaks into domains 
of different orientation (Seynaeve et al. [252]). The exact transition from single domain to multi-domain 
structures depends on the shape, dimensions, temperature and particular material, among other parameters. 
More recently, Villegas et al. [2151 1253j and Hoffmann et al. |254j experimentally investigated the switching 
of the ferromagetic dots from single domain to magnetic vortex state while sweeping the external field and the 
influence of their stray fields on the resistivity of the S/F hybrid sample (Fig. [^5)1 . The interaction between a 
vortex in a superconducting film and a magnetic nanodisk in the magnetic vortex state was studied theoretically 
by Carneiro [255] . For magnetic dots big enough to host a multi-domain state it is possible to tune the average 
magnetic moment by partially magnetizing the sample in a field lower than the saturation field or even recover 
the virgin state by performing a careful degaussing procedure similar to that shown in Fig. [5] (Gillijns et al. 
[TOl 1123] . Lange et al. |256| ). Interestingly, recently Cowburn et al. [178J showed that small disks of radius 
about 50 nm made of supermalloy (Ni-80%, Fe-14%, Mo-5%) lie in a single domain state with the magnetization 
parallel to the disk plane and with the property that their direction can be reoriented by small applied fields. This 
systems represents the closest experimental realization of in-plane free-rotating dipoles, which was theoretically 
analyzed by Carneiro |162l 11631 1164] within the London formalism. 

Whatever the mechanism of pinning produced by the magnetic dots, either core or electromagnetic, it 
is now a clearly established fact that changing the domain distribution in each dot has profound effects on 
the superconducting pinning properties as demonstrated, for example, by van Bacl et al. |196j . van Look et 
al. [202] . This result points out the importance of performing a careful study of the magnetic properties of 
the dots in order to identify the domain size, shape, distribution, and stable states. Van Bael et al. |196j 
presented the first report directly linking changes in the hysteresis loop of a superconducting Pb film when 
the underlying submicron Co islands are switched from 2x2 checkerboard magnetic domain pattern to single 
domain structures. 

As we pointed out above, unfortunately both, the multidomain state and the magnetic vortex state, 
still involve a sizable component of the magnetic stray field which eventually influences the response of the 
superconducting properties by locally suppressing the order parameter. In other words, it is actually not 
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Figure 26. (color online) Different magnetic states realized in square permalloy micro-loops depending on the 
direction of the applied magnetic field, by courtesy of Metlushko et al. 



possible to completely switch off the magnetic pinning using singly connected structures. It has been recently 
shown that a way to partially circumvent this limitation can be achieved by using multiply connected ring-like 
magnetic structures. In this case, if a flux-closure state is induced in the magnetic ring, in principle, there is 
almost no stray field present, besides small fields due to domain walls at the sharp corners of the ring. Indeed, 
a two-dimensional magnetic material of ring like shape of group symmetry C„ can be set in two flux closure 
states of opposite chirality and n{n — 1) different polarized states. In a square loop, for instance, 12 states 
corresponding to six different dipole directions with two opposite dipole orientations are expected. If the net 
dipolar moment is parallel to the side of the square, the final state is named horseshoe state whereas if the 
dipole is along the diagonal of the square, it is called onion state. Figure [25] shows the different topologically 
non-equivalent magnetic states for a square ring of magnetic material with in-plane magnetization. Clear 
experimental evidence of ON/OFF magnetic pinning potentials induced by these type of multiply-connected 
ferromagnetic structures have been demonstrated by Silhanek et al. \257\ \258\ I259j . It is worth mentioning that 
the S/F structures investigated in Ref. |258j exhibit two well distinguished phases corresponding to a disordered 
phase when the sample is in the as-grown state and an ordered phase when the sample is magnetized with an 
in-plane field. These order-disorder transitions manifest themselves as an enhancement of submatching features 
in the field dependence of the critical current which cannot be explained from a simple rescaling of the response 
corresponding to the disordered phase. 

Random (disordered) magnetic inclusions 

Early studies of the influence of ferromagnet on the superconducting state were performed in the sixties 
by Strongin et al. |260j . Alden and Livingston [2611 1262j and Koch and Love [263| for a dispersion of fine 
ferromagnetic particles (Fe, Gd, Y) in a superconducting matrix. These reports motivated further experimental 
and theoretical investigations of the influence of randomly distributed particles on/underneath/inside 
superconducting materials, which continue nowadays (Sikora and Makiej [264|, [265], Wang et al. [266], 
Lyuksyutov and Naugle [2671 EMI |26g, Santos et al. [2TO], Kuroda et al. [271], Togoulev et al. [272] . 
Kruchinin et al. [273], Palau et al. |274l [275], Haindl et al. [2761 [277], Snezhko et al. [278], Rizzo et al. 
[HI], Stamopoulos et al. [IHUl [Ml [Ml [IHUM] , Suleimanov et al. [MS], Xing et al. [Mi[M3IMS])- In 
most of these investigations no precautions were taken to electrically isolate the magnetic particles from the 
superconducting material which presumably results in a substantial core pinning due to proximity effects. 

Xing et al. [287j reported on controlled switching between paramagnetic0 and diamagnetic Meissner effect 
in S/F nanocomposites consisting of Pb films with embedded single-domain Co particles. These authors argue 
that in this particular system the paramagnetic Meissner effect attributed to the superconducting part only, 
originates from the spontaneous formation of vortices induced by the ferromagnetic inclusions. Therefore, the 
different contributions of the external field and the spontaneous vortices to the resulting magnetization of the 
sample, make it possible to manipulate the sign of Meissner effect by changing the orientation of the magnetic 



t The paramagnetic Meissner effect in various superconducting systems is discussed in the review of Li 12891 . 
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Figure 27. (color online) Bottom parts of the magnetization curves M vs. Hext for a superconducting Pb film 
covering a Co/Pt multilayer. The curves M{Hext) corresponding to the different values of the parameters s, 
which is defined as the fraction of magnetic moments that are pointing up (m > 0) relative to the total number 
of magnetic moments: s = 0.1 (open circles), s = 0.3 (diamonds), s = 0.5 (crosses), s = 0.7 (triangles) and 
s = 0.85 (filled circles), adapted from Lange et al. |256| . 



moments embedded in the supercondicting matrixlj. 

Vortex dynamics in a periodic magnetic field 

Here we want to briefly discuss the pecuharities of low-frequency vortex dynamics in nonuniform magnetic fields. 
Magnetic templates placed in the vicinity of a superconducting film not only induces changes in the static pinning 
properties but also in the overall dynamic response of the system. Lange et al. |206j demonstrated that the 
vortex-antivortex pairs induced by an array of out-of-plane magnetized dots lead to a strong field polarity 
dependent vortex creep as evidenced in the current- voltage characteristics. This result shows that in S/F 
hybrids with perpendicular magnetized dots vortices and antivortices experience a different pinning strength. 
A theoretical study of the dynamic evolution of these interleaved lattices of vortices and antivortices in the case 
of in-plane point like dipoles has been recently addressed by Carneiro [291] and Lima and de Souza Silva [292j . 

A more subtle effect, namely magnetic-dipole-induced voltage rectification was predicted by Carneiro [243j . 
Unlike conventional ratchet systems|§|, in the particular case of magnetic ratchet, induced by in- plane magnetized 
dots, the motion of vortices is in the opposite direction than the motion of antivortices, thus giving rise to a 
field-polarity independent rectification (Silhanek et al. [2581 1259] , de Souza Silva et al. |303j V This magnetic 
dipole-induced ratchet motion depends on the mutual orientation and strength of the local magnetic moments 
thus allowing one to control the direction of the vortex drift. In some cases, a nonzero rectified signal is observed 
even at Hext = resulting from the interaction between the induced vortex-antivortex pairs by the magnetic 
dipoles [303] . It is worth emphasizing that in the case of in-plane magnetic dipoles treated by Carneiro '243j. 
the inversion symmetry is broken by the stray field of the dipoles, thus giving rise to different depinning forces 
parallel and anti-parallel to the dipoles orientations, as shown in Fig. 

I Previously, Monton et al. 12901 reported on an experimental observation of the paramagnetic Meissner effect in Nb/Co 
superlattices in field-cooled measurements, however the origin of this effect remains unclarificd. 

§ Early theoretical studies showed that a vortex lattice submitted to an oscillatory excitation in the presence of a non-ccntro- 
symmetric pinning potential gives rise to a net drift (v) of the vortex lattice which in turn generates a dc voltage signal 
Vdc = _/ [(v) X Hea;i] ■ dl along the direction of bias current (Zapata et al. 12931 . Lee et al. 12941 . Wambaugh et al. 12951 ). 
These predictions are in agreement with recent experimental results obtained for purely superconducting systems (Villegas et al. 
[296], Wordenweber et al. [297] , van de Vondel et al. [2981 . Togawa et al. [2991 . de Souza Silva et al. [300) . Wu et al. [30T] . 
Aladyshkin et al. [302]'). 
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3.4. Planar S/F bilayer hybrids 

In this section we shall discuss the properties of continuous planar S/F structures which have macroscopically 
large lateral dimensions. As before, the superconducting and ferromagnetic films are assumed to be electrically 
insulated from each other. 

Appearance of vortices in planar S/F structures 

The interaction of the Meissner currents and the currents induced by vortex lines with a one-dimensional 
distribution of the magnetization (both single domain walls, periodic domain structures and magnetic bars) in 
the London approximation was considered by Sonin [53] , Genkin et al. [304] , Bespyatykh and Wasilevski [305j , 
Bespyatykh et al. [306], Helseth et al. [307], Laiho et al. [308], Traito et al. [309], Erdin [310|, Bulaevskh 
and Chudnovsky [3111 1312] , Kayali and Pokrovsky [313) , Burmistrov and Chtchelkatchev [314j , Ainbinder and 
Maksimov [315j . Maksimova et al. [3161 1317] . It was found that in order to create vortex-antivortex pairs 
in the S/F bilayer with out-of-plane magnetization at Hf,r^t = (and thus keeping the total flux through the 
superconducting film zero) the amplitude of the magnetization AIs should overcome the following threshold 
value [3051130811309] 

My-av = -j oc — — InA/C, 

4a w X^w 

where a is a numerical factor of the order of unity and the period w of the domain structure is assumed 
to be fixed (the hard-magnet approximation). This estimate corresponds to the case when the width of the 
domain walls is much smaller than other relevant length scales. The critical magnetization M^_^^ decreases 
monotonically with decreasing superconducting film thickness Dg. The equilibrium vortex pattern appearing in 
the superconducting film at Hext = and Mg > M^^y consists of straight vortices, arranged in one-dimensional 
chains, with alternating vorticities corresponding to the direction of the magnetization in the ferromagnetic 
domains [3051 13081 1309| . The parameters of such vortex configuration with one or two vortex chains per half- 
period was analyzed by Erdin [310] , It was shown that in equilibrium the vortices in the neighboring domains 
are halfway shifted, while they are next to each other in the same domain. Alternatively, as the thickness Dg 
increases, the vortex configuration, consisting of vortex semi-loops between the ferromagnetic domains with 
opposite directions of the magnetization, becomes energetically favorable [30811309] provided that Mg > M^j^^p^, 
where 

Hci $0 InA/e 

JW7„_„„ — — ; ; — CX 



'°°P' 8\n{w/nX) A2 ln(w;/^A)' 
The destruction of the Meissner state in the S/F bilayer with in-plane magnetization was considered by 
Burmistrov and Chtchelkatchev [314] . Since the out-of-plane component of the field, which is responsible for 
the generation of the vortex, is maximal near the domain wall (unlike from the previous case) and goes to zero 
in the center of magnetic domains, one can consider only a single domain wall. At Hext = a creation of a 
single vortex near the Bloch-type domain wall of width d corresponds to the condition 
Hci A r 2A/(5, 7r5/(4A) < 1 



X 



47r Df il-32A/(7r2,5), 7r(5/(4A) > 1. 



Magnetic pinning and guidance of vortices in planar S/F structures 

Irrespective of whether the domain structure in the ferromagnetic layer is spontaneously created or was present 
beforehand, the spatial variation of the magnetization will lead to an effective vortex pinning (Bespyatykh et 
al. [306] . Bulaevskii et al. [318] V However, there are discrepancies in the estimates concerning the pinning 
effectiveness. Indeed, Bulaevskii et al. [318[ argued that superconductor/ferromagnet multilayers of nanoscale 
period can exhibit strong pinning of vortices by the magnetic domain structure in magnetic fields below the 
coercive field when the ferromagnetic layers exhibit strong perpendicular magnetic anisotropy. The estimated 
maximum magnetic pinning energy for a single vortex in such a system is about 100 times larger than the 
core pinning energy produced by columnar defects. In contrast to that, Bespyatykh et al. L306| have shown 
that the effectiveness of magnetic pinning of vortices in a layered system formed by an uniaxial ferromagnet, 
does not considerably exceed the energy of artificial pinning by a column- type defect, regardless the saturation 
magnetization of the ferromagnet. The limitation of the pinning energy is caused by the interaction of external 
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vortices with the spontaneous vortex lattice formed in the superconducting film when the magnetization of the 
ferromagnetic film exceeds the critical value (see Eq. 

There have been numerous experimental investigations corroborating the enhancement of the critical current 
in planar S/F hybrids. It was shown that the presence of a bubble domain structure in Co/Pt ferromagnetic 
films with out-of-plane magnetization modifies the vortex pinning in superconducting Pb films (Lange et al. 
[2051 12561 13191 [52n] l. leading to an increase of the width of the magnetization loop AI{Hext) as compared with 
a uniformly magnetized S/F sample (Fig. I27p . The crossover between an enhanced magnetic pinning on bubble 
magnetic domains observed at low temperatures and a suppressed magnetic pinning at temperatures close to 
Tc for a demagnetized S /F bilayer can be possibly associated with an increase of an effective penetration length 
/Ds, characterizing the vortex size, and an effective averaging on the small-scale variation of the nonuniform 
magnetic field provided that X^/Dg considerably exceed the period of the magnetic field (Lange et al. [320] 1. 
Interestingly, the parameters of the bubble domain structure (the size and the density of domains of both signs 
of magnetization) can be controlled by demagnetization similar to that reported in Refs. [96lll23j . A three-fold 
enhancement of the critical depinning current in Nb films fabricated on top of ferromagnetic Co/Pt multilayers 
was observed by Cieplak et al. |3211 1322] based on magnetization measurements and on the analysis of the 
magnetic field distribution obtained by using a ID array of Hall sensors. The mentioned enhancement of the 
magnetic pinning takes place in the final stages of the magnetization reversal process, and it can attributed to 
residual un-inverted dendrite-shaped magnetic domains. 

High-resolution magneto-optical imaging performed by Goa et al. |323j in superconducting NbSe2 single 
crystals and ferrite-garnet films demonstrates that the stray field of Bloch domain walls can be used to 
manipulate vortices. Indeed, depending on the thickness of the sample, the vortices are either swept away 
or merely bent by the Bloch wall. 

Vlasko-Vlasov et al. |247[ 1324] and Belkin et al. j325[ 1326] studied the anisotropic transport properties of 
superconducting MoGe and Pb films and NbSe2 single crystals which are in the vicinity of a ferromagnetic 
permalloy film. In these works a quasi-one-dimensional distribution of magnetization can be achieved by 
applying a strong enough in-plane field H^xt > 300 Oe, which aligns the domain walls in a desired direction. 
Such a domain structure was maintained even after switching off the external magnetic field. Magneto-optical 
measurements directly display the preferential direction of the vortex entry in the presence of the perpendicular 
magnetic field along the domain walls [the panel (c) and (d) in Fig. [23]. By reorienting the magnetic domains 
using combinations of dc and ac fields, it is possible to rearrange current patterns in the S/F bilayer and thus 
manipulate its conductivity. The presence of this rotatable periodic stripe-like magnetic domain structure 
with alternating out-of-plane component of magnetization results in a difference in the critical depinning 
current density between cases when the magnetic domain stripes are oriented parallel and perpendicular to 
the superconducting current: j| > J^. For planar thin-film Pb/Py structures Vlasko-Vlasov et al. |324j 
observed a pronounced magnetoresistance effect yielding four orders of magnitude resistivity change in a few 
millitesla in-plane field. In addition, the S/F bilayer exhibits commensurability features that are related to the 
matching of the Abrikosov vortex lattice and the magnetic stripe domains (Belkin et al. |326| ). The matching 
effects are less apparent than for S/F hybrids with magnetic dots, although commensurability becomes more 
pronounced as temperature is lowered. This result can be explained by the gradual decrease in the A value, 
which leads to stronger modulation of the magnetic field in the superconductor at lower temperatures and 
consequently, to more prominent magnetic interaction with ferromagnetic domain structure. 

It is interesting to note that the effect of magnetic domains on the pinning of vortices, was also observed 
in high— Tc superconductors such as YBa2Cu307_5 (Garcfa-Santiago et al. [327], Jan et al. [328], Zhang et 
al. [329] ■ Laviano et al. [330j ). At the same time the inffuence of the ferromagnet on the nucleation in the 
high— Tc superconductors should be rather small due to extremely small coherence length (of the order of few 
nanometers) . 

Current compensation effect and field-polarity dependent critical current 

A superconducting square with in-plane magnetized ferromagnet on top was proposed by Milosevic 
et al. |331j as a potential field and current compensator, allowing to improve the critical parameters of 
superconductors. Indeed, such a magnet generates stray fields of the same amplitude but opposite signs at 
the poles of the magnet, therefore the field-compensation effect leads to the enhancement of the upper critical 
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Figure 28. (color online) Diode effect in the Nb bridge (2 /xm width) with the in-plane magnetized Co stripe 
on top: Experimental dependence of (the critical current in the x— direction, see the geometry of the S/F 
system on the inset) and 1^ (the critical current in the opposite direction) on the external magnetic field Hext 
applied in the y— direction (T = 4.2 K, Tco = 9-2 K), adapted from Vodolazov et al. 13331 . In the inset the dc 
I — V characteristic of our hybrid system H^xt = 500 Oc is presented, showing a pronounced diode effect. 



field equally for both polarities of the external field. The superconducting state was shown to resist much higher 
applied magnetic fields for both perpendicular polarities. In addition, such ferromagnet induces two opposite 
screening currents inside the superconducting film plane (in the perpendicular direction to its magnetization), 
which effectively compensates the bias current, and therefore superconductivity should persist up to higher 
applied currents and fields. These effects have been recently studied experimentally by Schildermans et al. 
[332j in an Al/Py hybrid disk of 1.7 /xm diameter where a finite dipolar moment lying in the plane of the 
structure was achieved by pinning magnetic domains with the contact leads used for electrical measurements. 

Vodolazov et al. [333] and Touitou et al. [334] considered an alternative experimental realization of the 
current compensator, consisting of a superconducting bridge and a ferromagnetic bar magnetized in-plane and 
perpendicularly to the direction of the bias current. Such geometry allows one to weaken the self- field of 
the superconducting bridge near its edge and thus to enhance the total critical current corresponding to the 
dissipation-free current flow. Since the self-field compensation occurs only for a certain direction of the current 
(for fixed magnetization), the presence of magnetized coating leads to a diode effect - the current-voltage 
I ^ V dependence becomes asymmetrical (Fig. [28]). Later the similar difference in critical currents flowing in 
opposite directions was studied experimentally by Morelle and Moshchalkov [335j for a system consisting of a 
superconducting Al strip placed close to a perpendicularly magnetized Co/Pd rectangle and Vodolazov et al. 
[336j for Nb/Co bilayer in the presence of titled external magnetic field. 

3.5. Stray field-induced Josephson junctions 

Josephson junctions consist of weak links between two superconducting reservoirs of paired electrons. 
Commonly, these junctions are predefined static tunnel barriers that, once constructed, can no longer be 
modified/tuned. In contrast to that, a new concept of the Josephson junctions with a weak link generated 
by the local depletion of the superconducting condensate by a "magnetic barrier" from a micro/nano-patterned 
ferromagnet can be realized (Sonin 98]). Interestingly, this type of devices offer an unprecedented degree of 
flexibility as it can be readily switched ON/OFF by simply changing between different magnetic states using 
an in-plane field. This switching process is fully reversible, and non-volatile since does not require energy to 
keep one of the magnetic states. 

A pioneer investigation of the properties of superconducting weak links achieved by local intense magnetic 
fields was performed by Dolan and Lukens [337]. The sample layout used by these authors and their typical 
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(b) Unquenched state 




Figure 29. (color online) (a) Sample layout investigated by Dolan and Lukens [337) : An uniform Al bridge 
was covered with a superconducting Pb strips (dashed rectangles) everywhere but in a small region near the 
center of the bridge. Due to the flux expulsion from the Pb strips the local magnetic field is primarily confined 
to this gap. 

(b) Sample configuration investigated by Clinton and Johnson 13381 l339l : a Pb (or Sn) transport bridge is 
partially covered with a ferromagnetic Py strip with in-plane magnetic moment M. When M is parallel to the 
bridge a strong stray field depletes the superconducting order parameter in a small region near the border of 
the Py bar (quenched state) thus inducing a weak link. 

(c) Schematic presentation of magnetoquenched superconducting quantum interferometer, consisting of two 
superconducting Pb bridges connected in parallel and permalloy film on top, after Eom and Johnson |342| . 

dimensions is schematically shown in Fig. [551 (a) and it consists of an Al bridge locally covered by a plain Pb 
film which has a thin gap of width g and spans the width of the Al strip at its center. By applying an external 
magnetic field, the Pb film screens the magnetic field due to the Meissner effect in the whole Al bridge but 
magnify its intensity at the gap position. This effect leads to a local region of suppressed superconductivity 
which gives rise to dc and ac Josephson effect as evidenced by a finite critical current and the presence of Shapiro 
steps in the current-voltage characteristics at = nfioj / (2e) when the system was irradiated with rf-excitations 
with frequency lo = 2Trf, n is integer. Interestingly, the Josephson-like features appear for applied fields in the 
shield gap approximately equal to the upper critical field of the Al film. 

An alternative method to obtain a field-induced weak link has been more recently introduced by Clinton 
and Johnson |338[ 13391 1340[ I341j . The basic device consist of a bilayer of a thin superconducting strip and a 
ferromagnetic layer with in-plane magnetic moment overlapping the width of the bridge [see panel (b) in Figure 
129] . When the magnetic moment is parallel to the superconducting bridge the dipolar fringe is strong enough to 
locally suppress the superconducting order parameter across the bridge (quenched state) and thus create a weak 
link. This effect can be turned off by simply magnetizing the ferromagnetic layer perpendicular to the transport 
bridge with an external in-plane dc field or by a current pulse in a separate transport line [341] . Clearly, the 
proposed switchable Josephson junction seems to be very attractive for potential technological applications, 
since energy is required only to change the magnetic states, which are thereafter maintained in thermodynamic 
equilibrium. Later on based on the same idea Eom and Johnson [342] proposed a switchable superconducting 
c^uantum interferometer consisting of a ferromagnetic Py film partly covering two parallel superconducting Pb 
bridges fabricated in a loop geometry. The dependence of the voltage V , induced on this superconducting loop 
at injection of stationary bias current, on the perpendicular magnetic field Hext is shown in Fig. I30f b) and it 
reminds the standard Fraunhofer diffraction pattern (Barone and Paterno [34j ) . 
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Figure 30. (color online) (a) The I — V curves obtained for a plain superconducting Pb bridge (2 fj,m wide) 
subjected in the inhomogeneous magnetic field, quenched state [see the panel (b) in Fig. I29j for different 
intensities of rf-irradiation, adapted from Clinton and Johnson 13391 . The experiment was carried out at T = 5 
K, T/TcO ~ 0.76, Hext = 0, frequency / of the radio signal equal to 0.75 GHz. 

(b) The I — V dependence obtained for a superconducting Pb loop of the width 4.5 fim with rectangular hole 
1.5x7.0 fircP covered by permalloy film [see the panel (c) in Fig. 1291 . adapted from Eom and Johnson |342l . This 
curve demonstrates the short period oscillations with the period determined by the area of the superconducting 

loop Sloop- 
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4. Hybrid structures: superconductor soft magnets 

Thus far, we have discussed tlic influence that a ferromagnet has on the superconducting properties of S/F 
hybrids, assuming that the magnetization of the ferromagnet M remains practically unaltered. In this last 
section, we consider the possibility that the magnetization M can be changed either due to the external magnetic 
field or by the superconducting screening currents induced by the magnetic subsystem, which are particularly 
relevant at low temperatures. This situation could, in principle, be achieved by using paramagnetic materials 
or soft ferromagnetic materials with a low coercive field. 

The equilibrium properties of "superconductor - soft magnet" hybrid structures (so-called soft S /F hybrids) 
can be obtained phenomenologically by the minimization of the Ginzburg-Landau energy functional Eq. ([l]) or 
the London energy functional Eq. ([T7]), in which the term G,„ responsible for the self-energy of the ferromagnet 
becomes important 

G,n = y (j|| |Vil4|2 + J,i |V7\4|2 + iVAf.p) dV- j 2nQMl dV, (23) 

" Vf Vf 

where J\\ and J_l characterize the exchange interaction between spins in a uniaxial ferromagnet with respect to 
the in-plane and out-of-plane direction, Q is a quality factor taking into account the internal anisotropy of the 
ferromagnet and determining the preferable orientation of the magnetization (either in-plane or out-of-plane). 
Equation (j23p describes the energy cost for having a slowly varying spatial distribution of the magnetizatior|jJ 
and, in particular, it describes the energy of a domain wall in a ferromagnet. In some cases (for instance, for 
rapid M variations typical for ferromagnets with domain walls of rather small width), in order to simplify the 
problem, the increase of the free energy given by Eq. ((23|) can be taken into account phenomenologically by 
substituting Gm by a fixed term Gdw representing the energy of a domain wall. 

Modification of the domain structure in a ferromagnetic film by the superconducting screening currents 
The influence of superconducting environment (both substrate or coating) on the equilibrium width of magnetic 
domains in ferromagnetic films was considered theoretically by Sonin [Ml, Genkin et al. [304] . Sadreev [344j . 
Bespyatykh et al. [345[ I346j , Stankiewicz et al. |347[ 1348] , Bulaevskii and Chudnovsky [311[ 1312] , Daumens 
and Ezzahri [349]. In particular, one can expect a prominent change in the equilibrium period of a one- 
dimensional domain structure at H^xt — for rather thick ferromagnetic films {Df 3> w) with out-of-plane 
magnetization. Indeed, the Meissner currents, induced by the ferromagnet, will decrease the magnetic field 
inside the superconductor (usual fiux expulsion effects) and significantly increase the magnetic field inside the 
ferromagnet. As a consequence, the density of the free energy of the ferromagnet, given by B^/Stt — B • M or, 
equivalently by H^/Stt - 2ttM'^, raises for a given Mz distribution. However, the total energy of the S/F system 
can be lowered by a decrease of the period of the ferromagnetic domains: the smaller the period, the faster the 
decay of H away from the surfaces of the ferromagnetic film. Thus, it is expected that the equilibrium width of 
magnetic domains in planar S /F bilayer becomes smaller below the critical temperature of the superconducting 
transition as compared with the state T > Tcq. In contrast to that, for thin ferromagnetic films {Df ^ w) the 
opposite behavior is predicted: the domain width in the free ferromagnetic film should be smaller than that for 
the same film on top of a superconducting substrate [348] . This can be understood by taking into account the 
change of the far-zone demagnetizing field characteristics. In addition, Stankiewicz et al. |348j argued that the 
effect of the superconducting substrate on the period of domain structure in ferromagnetic films with in-plane 
magnetization is rather small as compared with that for the out-of-plane magnetized ferromagnets. However, 
an increase of the magnetostatic energy of the S /F hybrids at T < Tco due to the Meissner currents results also 
in a shrinkage of the equilibrium width of an isolated 180° Bloch wall, separating two ferromagnetic domains 
with in-plane magnetization, in the vicinity of the superconducting substrate, as was predicted by Helseth et 
al. [307] . 

The foreseen decrease of the period of the domain structure in a ferromagnetic garnet film in contact with 
a superconducting Pb film was recently investigated by magneto-optical imaging (Tamegai et al. (350] ). It was 
demonstrated that the shrinkage depends both on temperature and the thickness of superconducting coating 

f The theory of superconductor - soft ferromagnet systems near the "ferromagnet-paramagnet" transition has been considered by 
Li et al. |343) within Ginzburg-Landau formalism. 
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Pb film on top magnetic garnet film 




Figure 31. (color online) (a) Tlie temperature dependence of tiie slirinkage ratio of tlie width of magnetic 
domains in a hybrid structure consiting of a Pb film {Ds = 300 nm) on top of a magnetic garnet film in 
comparison with the same ferromagnetic film without superconducting coating, adapted from Tamegai et al. 

(b) The magneto-optical image of four segments of the superconducting Pb/garnet film structure, differing by 
the thickness of the Pb film (D^ = nm, 100 nm, 200 nm, 300 nm in a counter-clock- wise direction), taken at 
T = 5.0 K, adapted from Tamegai et al. |350| . 

layer. The temperature dependence of the shrinkage factor s evaluated by comparing the average width (w) of 
the magnetic domain width in regions with and without the superconducting Pb film is shown in Fig. ISlT a). It 
points out that the lower the temperature, the narrower the magnetic domains are (s = 0.47 at T = 5.0 K). 

Alteration of magnetization of ferromagnetic dots by the superconducting screening currents 
The Meissner currents also influence the magnetic states and the process of magnetization reversal in 
ferromagnetic disks placed above a superconductor. It is well known that a uniformly magnetized (single 
domain) state is energetically favorable for radii Rf smaller than some threshold value Rj (for a given thickness 
of the dot Df), while the magnetic vortex state can be realized for Rf > Rj. The typical M{Hext) dependence 
for ferromagnetic disks for i?/ > Rj was already shown in Fig. [551 The dependence Rj vs. Df (the phase 
diagram in the "diameter-height" plane) in the presence of a bulk superconductor, characterized by the London 
penetration depth A, was investigated numerically by Fraerman et al. |351j and later analytically by Pokrovsky 
et al. [352] ■ It was shown that the smaller A, the smaller the critical diameter RJ becomes for a given dot's 
thickness. The transitions between the two magnetic states can be induced also by increasing the external 
magnetic field: the magnetic vortex, possessing an excess energy at zero field, becomes energetically favorable 
for finite external fields (the magnetic- vortex nucleation field). Although the energy of the interaction between 
the superconductor and ferromagnet is expected to be much smaller than the self-energy of the ferromagnetic 
particle (for realistic A values), it could lead to an experimentally observable decrease in the magnetic vortex 
nucleation field ^^'^ increase in the magnetic vortex annihilation field Hann (Fig- [3H)- 

The appearance of a spontaneous magnetization of individual S/F hybrids, consisting of an Al bridge and 
demagnetized Ni dots on top, upon cooling through the superconducting transition temperature at Hext = 0, was 
reported by Dubonos et al. [353]. Indeed, the reshuffling of magnetic domains in the submicron ferromagnetic 
disk, caused by temperature-dependent screening of the domain's stray fields by the superconductor, can explain 
the observed appearance of nonzero magnetization of the ferromagnet at low temperatures. More recently, the 
modification of the magnetic state of Nb/Co and Nb/Py superlattices induced by screening currents in the 
superconducting Nb films was studied experimentally by Monton et al. |2901 13541 1355] and Wu et al. [356] . 

Kruchinin et al. [273] demonstrated theoretically that a superconducting environment modifies the 
magnetostatic interaction between localized magnetic moments (embedded small ferromagnetic particles), 
resulting either in parallel or antiparallel alignment of neighbor dipolar moments at -ffgait = 0- The crossover 
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between these regimes depends on the ratio of the interparticle spacing and the London penetration depth, 
and thus preferable "magnetic" ordering (ferromagnetic vs. antiferromagnetic arrangements) can be tuned by 
varying temperature. 

Mixed state of soft S/F hybrid structures 

The magnetostatic interaction between a vortex-free superconducting film and a uniformly magnetized 
ferromagnetic film at Hext — may cause the spontaneous formation of vortices in the superconductor 
and magnetic domains in the ferromagnet in the ground state of planar S/F bilayers with perpendicular 
magnetization. Lyuksyutov and Pokrovsky [357j . and Erdin et al. |358j argued that the ground state of 
the S/F system could be unstable with respect to the formation of superconducting vortices. Indeed, for a 
uniformly magnetized S/F bilayer, characterized by a magnetization of the ferromagnetic film per unit area 
m — MgDf, the magnetostatic interaction between the superconductor and the ferromagnet changes the total 
energy of an isolated vortex line to = ei"' — m$o |146j as compared with the self-energy of the vortex in the 
superconducting film without a ferromagnetic layer. As a consequence, the formation of vortices becomes 
energetically favorable as soon as e^, < (either for rather large m values or at temperatures close to Tcq where 
el,"' vanishes). However, as the lateral size of the S/F system increases, the averaged vortex density riy would 
generate a constant magnetic field Bz — riy^Q along the z— direction which can lead to an energy increase 
larger than the gain in energy due to creation of vortices. Hence, in order for the vortex phase to survive, 
the ferromagnetic film should split in domains with alternating magnetization in a finite temperature range at 
T < TcQ. As long as the magnetic domain width exceeds the effective penetration depth, the energy of the 
stripe domain structure seems to be minimal (Fig. I33|) . Interestingly, the interaction between a single vortex 
in a superconducting film and the magnetization induced by this vortex in the adjacent ferromagnetic film can 
cross over from attractive to repulsive at short distances (Helseth |171j V 

Carneiro studied the interaction between superconducting vortices and a superparamagnetic particle with 
constant dipolar moment, which is assumed to be able to freely rotate, in the London model [1621 1163"! I164j . 
It was found that, due to the rotational degree of freedom, the pinning potential for superconducting vortices 
differs significantly from that for a permanent dipole. In particular, the interaction between the superconducting 




Figure 32. (color online) (a) The dependence of the magnetic vortex annihilation field Ha„„, corresponding to 
the transition from the magnetic vortex state to a single-domain state, on the diameter of disk 2i?y, calculated 
for an isolated ferromagnetic disk of 20 nm thickness (i.e, without superconductor, A = oo, open circles) and 
for the same disk placed above a bulk superconductor (A = 50 nm, filled circles), adapted from Fraerman et al. 
[351]. 

(b) The dependence of the magnetic vortex nucleation field ffJI^^,;, corresponding to the transition from single- 
domain state to the magnetic vortex state for the same problem, adapted from Fraerman et al. 351 . 

(c) The magnetization curve AfH/Ms vs. in-plane external field -ffj^j demonstrating the process of the 
magnetization reversal for the magnetic disk (20 nm thickness and 100 nm diameter) for A = oo (open circles) 
and A = 50 nm (filled circles), adapted from Fraerman et al. |351) . Thus, the screening effect increases the 
width of the Hext interval at the ascending and descending branches of the magnetization curve where magnetic 
vortex state is energetically favorable. 
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vortex and the magnetic dipole can be tuned by applying an in-plane external field: the corresponding depinning 
critical current was shown to be anisotropic and its amplitude potentially varies by as much as one order of 
magnitude. Later on, this approach was generalized by Carneiro [255j for hybrid systems consisting of thin 
superconducting film and a soft ferromagnetic disks in the magnetic vortex state (similar to that for Refs. 
[3511 1352] but considering a vortex line inside the superconducting sample) . 

A new method of pinning vortices in S/F epitaxial composite hybrids consisting of randomly distributed 
Gd particles incorporated in a Nb matrix was reported by Palau et al. [2751 1274] . Since the size of Gd particles 
are much smaller than the coherence length and the interparticle distance is much shorter than the penetration 
depth, this regime of collective magnetic pinning differs both from conventional core and magnetic pinning 
mechanisms. In this case, since a vortex "feels" a homogeneous superconductor (for length scales on the order 
of A), pinning effects are expected to be small. However, due to the local field of a vortex, the Gd particles 
can be magnetized and a moving vortex would lead to hysteretic losses in the magnetic particles, which in turn 
results in an increased pinning (for decreasing magnetic fields). 

Superconductor - paramagnet hybrid structures 

An alternative way of modifying the superconducting properties of soft hybrid structures is by using 
paramagnetic constituents, characterized by zero or very low remanent magnetization. Such superconductor- 
paramagnet hybrids with a magnetization M = (/i — 1) H/47r depending on the external field (/i is the magnetic 
permeability) in the presence of transport current were considered theoretically by Genenko [359] . Genenko 
and Snezhko [360 , Genenko et al. |3611 13621 1363j for /z ^ 1. It was predicted that the paramagnetic material 
placed near superconducting stripes and slabs can drastically modify the current distribution in such hybrids, 
thus, suppressing the current enhancement near the superconducting sample's edges inherent for any thin-film 
superconductor in the fiux-free current-carrying state. As a consequence, the current redistribution leads to an 
increase of the threshold value of the total bias current corresponding to the destruction of the Meissner state. 
In other words, the magnetically shielded superconductors even in the Meissner state are able to carry without 
dissipation rather high transport current comparable with the typical current values for a regime of strong 
flux pinning [3591 13611 13621 1363j . A survival of the Meissner state for thin-film superconducting rings carrying 
a current and placed between two coaxial cylindrical soft magnets was studied by Genenko et al. [3641 1365] 
and Yampolskii et al. [366| . The similar problem concerning with the distribution of magnetic field inside and 
outside a superconducting filament sheathed by a magnetic layer, as well as the magnetization of such a structure 
in the region of reversible magnetic behavior in the Meissner state was considered by Genenko et al. |367| . The 
formation of the mixed state in various superconductor/paramagnet structures in the presence of transport 




Figure 33. (color online) (a) Uniform state of the S/F bilayer above the superconducting critical temperature, 
(b) Spontaneously formed magnetic domain structure and coupled chains of superconducting vortices with 
alternating vorticity at T < Tco, adapted from Erdin et al. |358| . Solid arrows correspond to the magnetization 
of the ferromagnet, while black arrows schematically show the circulating superconducting currents. 
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current, or an external magnetic field H^^t or the field of hard magnetic dipoles, were analyzed by Genenko 
et al. |368j , Genenko and Rauh [369] , Yampolskii and Genenko [370j , and Yampolskii et al. [3711 13721 1373] 
in the framework of the London model. The Bean-Livingston barrier against the vortex entry in shielded 
superconducting filaments was shown to strongly depend on the parameters of the paramagnetic coating and, 
as a result, the critical field at which the first vortex enters can be enhanced [368] . 

Since for the superconductor/paramegnet hybrids there are no changes neither in the vortex structure 
in the Meissner state of superconductor or in the magnetic state of paramagnetic elements, characteristics 
of superconductor-paramagnet hybrids are presumed to be reversible, and ac losses should be minimal for 
such structures. It stimulated the implementation of paramagnetic and ferromagnetic coatings in high- Tc 
superconductors in order to improve the critical current and reduce the ac-losses (Majoros et al. j374j , Glowacki 
et al. ^375j, Horvat et al. [37S[377], Duckworth ,378^, Kovac et al. [37S], Touitou et al. [331], Pan et al. [350] . 
Jooss et al. [MI], Gu et al. [552], Gomory et al. [5551 

Although, strictly speaking, permalloy is a ferromagnet with rather low-coercive field, it can behaves 
qualitatively similar to paramagnetic materials. Indeed, the magnetization vector for thin-film structures 
deviates from in-plane orientation if a perpendicular external field is applied. Such rotation of the magnetization 
of the dot toward the out-of-plane direction while sweeping the external field gives rise to a z— component of 
magnetization depending on the external field. The effect of the stray field generated by soft permalloy dots on 
the critical current Ic of the superconducting Al loops was experimentally studied by Golubovic and Moshchalkov 
[385j . The monotonous decrease in the period of the oscillation on the Ic{Hext) with increasing the Hext value 
was interpreted as a flux enhancement due to the increase of the out-of-plane component of the dot's magnetic 
moment. In this sense soft magnetic materials are promising candidates for the design of a linear magnetic flux 
amplifier for applications in superconducting quantum interference devices. 
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5. Conclusion 

We would like to conclude this work by formulating what, we believe, are the most exciting unsolved issues and 
possible interesting directions for further studies of S/F hybrid systems. 

Spontaneous formation of a vortex lattice and domain patterns. As we discussed in Section 21 the 
magnetostatic interaction between a vortex-free superconducting film and a uniformly and perpendicularly 
magnetized ferromagnetic film at zero external field, may lead to a spontaneous formation of vortex-antivortex 
pairs in the superconductor accompanied by alternating magnetic domains in the fcrromagnct in the ground 
state (see Fig. l33land Refs. [33ll358j ). To the best of our knowledge, thus far there are no experimental results 
confirming this prediction. In part, this is likely due to the difficulties of combining a very low coercive field 
ferromagnetic materials, needed to guarantee the free accommodation of magnetic domains, and a well defined 
out-of-plane magnetic moment. 

Magnetic pinning. Based on London equations we have clearly defined the magnetic pinning energy as 
mo • B„ (see section Since London equations are valid as long as core contributions are negligible, in 
principle this simple relationship holds for materials with k = A/^ ^ 1 and low temperatures. Unfortunately 
in the vast majority of the experimental reports so far, it seems that these conditions are not satisfied. Any 
other contributions such as local suppression of Tc, proximity effect, or local changes in the mean free path, 
which are not accounted for within London approximation, could lead to a deviation from the treatment in the 
framework of London model. The problem that remains unsolved so far is the identification of the most relevant 
mechanisms of vortex pinning in S/F hybrid systems. 

Thermodynamic properties of S/F heterostructures. Although the electric transport in superconduc- 
tor/ferromagnet hybrid systems has been intensively studied during the past decades, very little is known 
about their thermodynamic and thermal properties such as their entropy, specific heat, thermal conductivity, 
etc. From an academic point of view it would be very relevant to investigate the nature of the phase transitions, 
or present thermodynamic evidence of confinement of the superconducting order parameter. On the other hand, 
estimating the heat released when the system changes its state might also provide useful information for devices 
based on S/F heterostructures. 

Electromagnet- superconductor hybrids. Most of the research performed so far has been focussed on the 
effects of an inhomogeneous field generated by a ferromagnetic layer onto a superconducting film. As was early 
demonstrated by Pannetier et al. |97j in principle there is no difference whether this inhomogeneous field is 
the stray field emanating from a permanent magnet or the magnetic field generated by micro (nano) patterned 
current carrying wire on top of the surface of the film. The great advantage of the latter is the degree of 
flexibility and control in the design of the magnetic landscape and the external tuneability of its intensity. Such 
electromagnet-superconductor hybrids represent a promising alternative for further exploring the basic physics 
behind S/F hybrid. 

Time resolved vortex creation and annihilation. Josephson 7r-junctions give rise to spontaneous formation of 
half-integer flux quanta, so called semifluxons (Hilgenkamp et al. [386]). It has been theoretically demonstrated 
that for long Josephson junctions with zigzag 7r-discontinuity corners the ground state corresponds to a fiat phase 
state for short separation between corners whereas an array of semifluxons is expected for larger separations. 
Interestingly, by applying an external bias current it is possible to force the hopping of these semifluxons 
between neighboring discontinuities (Goldobin et al. [387j . This hopping of semifluxons could be identified 
through time resolve ac measurements with drive amplitudes above the depinning current. There are clear 
similarities between these arrays of semifluxons in zigzag Josephson systems and the vortex-antivortex arrays 
in S/F systems. Indeed, recently Lima and de Souza Silva [292] have shown theoretically that the dynamics of 
the vortex-antivortex matter is characterized by a series of creation and annihilation events which should be 
reflected in the time dependence of the electrical fleld. Experimental work corroborating these predictions are 
relevant for understanding the dynamics of flux annihilation and creation in S/F systems. 

All in all, the strive to comprehend the ultimate mechanisms ruling the interaction between ferromagnets 
and superconductors has made this particular topic an active theoretical and experimental line of research. 
We believe that these vigorous efforts will inspire further developments in this area of solid state physics and 
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Pb films 
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Table 1. Summary of experimental research on vortex matter in the S/F hybrids with dominant orbital 
interaction (suppressed proximity effect). 
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Table 2. Summary of theoretical research on vortex matter in the S/F hybrids with dominant orbital interaction 
(suppressed proximity effect) 
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